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Abstract: High energy cosmic neutrino fluxes can be produced inside relativistic jets 



Q^ \ oscillations of the mass states inside the Earth. We show that conversion has several new 

' features which are not realized in other objects, in particular interference effects ("L- and H- 

wiggles") induced by the adiabaticity violation. The v — v scattering inside jet and inelastic 
^ ' neutrino interactions in the envelope may produce some additional features at £J > 10^ 

■ GeV. We study dependence of the probabilities and flavor ratios in the matter-affected 

region on angles ^13 and 023) on the CP-phase 8, as well as on the initial flavor content 
and density profile of the star. We show that measurements of the energy dependence 
of the flavor ratios will, in principle, allow to determine independently the neutrino and 
astrophysical parameters. 
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^> ' under the envelopes of collapsing stars. In the energy range E ~ (0.3 — 10 ) GeV, flavor ' 

■ conversion of these neutrinos is modified by various matter effects inside the star and the CIZ) 

■ Earth. We present a comprehensive (both analytic and numerical) description of the flavor ^ 
. conversion of these neutrinos which includes: (i) oscillations inside jets, (ii) flavor-to-mass ^ — ^ 
, state transitions in an envelope, (iii) loss of coherence on the way to observer, and (iv) 
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1. Introduction 

It is difficult to overestimate importance of future detection of high energy cosmic neutrinos 
which will open new window to the Universe (for reviews see e.g. This detection will 
bring unique information on astrophysical sources of the neutrinos, as well as on neutrinos 
themselves: their propagation, interactions and flavor conversion. A number of different 
high-energy, > 1 GeV, neutrino sources have been proposed in literature, which includes 
active galactic nuclei (AGNs) |^, gamma ray bursts (GRBs) |^], core collapse supernovae 
(SNe) [Q], supernova remnants etc. Properties of neutrino fluxes, energy range, shape 
of the energy spectra and flavor content depend on physical conditions in the sources, as 
well as on effects of propagation between the source and the Earth. In particular, the flavor 
content of fluxes is modified by the vacuum oscillations between the production region and 
the Earth-based detectors. For known neutrino mass squared differences the oscillations 
are averaged due to large distances (baselines), and the oscillation effects are described by 
the averaged vacuum oscillation probabilities which do not depend on neutrino energy. 

The sources listed above are detected by their electromagnetic (EM) radiation with 
space- and ground-based detectors. No astrophysical source has been detected so far by its 
high energy neutrino signal. Apart from the sources visible in EM radiation, it is expected 
that there are various hidden sources of neutrinos: the sources in which the EM radiation 
produced in the inner part of the object is absorbed in the surrounding dense material and 
only neutrinos can pass through and be detected at the Earth. One such optically thick 
source is a core collapse SN with mildly relativistic jets emitted by a central engine, a black 
hole or a highly magnetized neutron star. These jets form promptly after the collapse ||6| 
and may not break the envelope unlike highly-relativistic GRB jets. Therefore they may 
not lead to emission of high energy non-thermal EM radiation. Observation of late time 
radio afterglow and explosion geometry provide evidences of a plausible hidden jet in some 
SNe 0. These hidden jets, however, can be the sources of high-energy neutrino fluxes 
generated by interactions of the shock accelerated protons with surrounding matter and 
radiation |, |8|,||. 

The number of hidden sources can be much larger than the number of observable ones, 
limited only by the ratio of type Ib/c and type II SNe to GRB rates. The combined SN 
rate from roughly 4000 galaxies within 20 Mpc is larger than one per year and the SN rate 
in nearby starburst galaxies such as M82 and NGC253 is much higher than in the Milky 
way 1^, Recent 7-ray observation of M82 and NGC253 with the Fermi Large Area 
Space Telescope also supports an increased SN activity . The prospect of high-energy 



neutrino detection from SNe in these galaxies by neutrino detectors such as IceCube[12| 



ANTARES H] and KM3NeT |lj is rather high, if a significant fraction of the SNe are 



endowed with mildly-relativistic jets. Optical follow-up triggered by neutrino events in 



IceCube along the direction of the neutrino trajectories will further enhance detection 
prospect by identifying the SN and extrapolating its light curve back to the explosion 
time m. 

Neutrinos produced in the hidden jets cross an envelope and in it was shown that 
in general matter effects on neutrino oscillations are not small. In particular, the minimal 
width condition (the lower bound on the column density of electrons on the way of neutrinos 
which is required for strong matter effect) |17] can be satisfied. It was shown that the matter 
effects inside a star can substantially modify the average vacuum oscillation probabilies in 
the range from 10^ GeV (the IceCube threshold) up to 10^ GeV. Neutrino fluxes at the 
surface of the Earth have been computed. As an observable in the forthcoming experiments, 
the ratio of shower events (induced mainly by the electron and tau neutrinos) and the muon 
track events (produced mainly by interactions of muon neutrinos and antineutrinos) has 
been proposed [l^. The ratio as a function of energy depends on values of the mixing angles 
^13) ^23) the CP-violation phase 6 as well as on density profiles of the stellar envelope. It was 
concluded, in particular, that with about 10^ events one can explore the neutrino properties 
such as the type of mass hierarchy and CP violating phase, provided that the mixing angles 
will be measured at future reactor and accelerator experiments with high accuracy. Also 
certain information about properties of the source can be extracted. The rate of detection 
of individual source from nearby galaxies is one in 5-10 years. If the neutrino telescopes 
are able to detect the diffuse flux of neutrinos from all the sources, still one can see some 
deviation of the flavor ratio from that produced by the vacuum oscillations. Inversely, 
observation of such a deviation in certain energy range can be explained by matter induced 
transformations and large population of the hidden sources. With high statistics this 
feature in the energy spectrum can be extracted from large background of atmospheric 
neutrinos. 

In this paper we have reconsidered the flavor conversion of neutrinos from hidden jets. 
Our results differ from those in Ref. \1G\- The difference originates from treatment of av- 
eraging and the coherence loss. Neutrinos produced by pions, muons and kaons in the 
strong magnetic field have very short wave packets and lose coherence due to separation 
of the wave packets So, neutrinos arrive at the surface of the Earth as incoherent 

fluxes of the mass eigenstates and coherence is not restored in detector. Therefore one 
should compute probabilities of the flavor-to-mass transitions (i/q — )• Ui) inside the star. In 
contrast, in Ref. [16| the flavor-to-flavor transitions have been computed from the produc- 
tion point to the surface of the star and then the oscillations on the way from the surface 
of the star to the earth have been averaged. This leads to different results of numerical 
computations. In particular, fast oscillatory behavior of probabilities with energy appears 
according to We have included in consideration also additional effects not considered 
in [ p^ ] , such as energy-dependent particle to antiparticle ratio and flavor ratio in the initial 
neutrino fluxes, oscillations in jets, inelastic interactions of neutrinos and neutrino- neutrino 
scattering. Detailed and comprehensive study of the conversion inside the star and on the 
way between the star and the Earth is performed. We show that the main effect is due to 
the adiabatic and partially adiabatic conversion (the MSW-effect) [^] inside an envelope. 
Both numerical and analytical results for probabilities, neutrino fluxes and flavor ratios at 



the Earth are presented. Computation of number of events in specific detectors is beyond 
the scope of this paper and will be given elsewhere. 

The paper is organized as follows. In Sec. 2 we describe the model of hidden source, 
summarize physical conditions at the neutrino production site, and present properties of 
the generated neutrino fluxes. We describe neutrino conversion inside the star in Sec. 3 
that includes the adiabatic conversion at low energies, adiabaticity violation at energies 
above the 1-3 resonance, interference effects in the range of adiabaticity violation which 
lead to H- and L- wiggles. We describe properties of the conversion probabilities in specific 
channels, their dependence on neutrino parameters and characteristic of density profiles in 
Sec. 4. In Sec. 5 we present neutrino fluxes and flavor ratios at the surface of the Earth. 
We consider their dependence on neutrino parameters, original flavor content as well as 
on the density profile of the stellar envelope. Conclusions are presented in Sec. 6. Some 
details of neutrino flux calculation, explanation of the difference of results of and this 
paper, and details of estimations of the Earth matter effects are presented in Appendices 
A, B and C correspondingly. 



2. Neutrinos from jets: production, physical conditions, fluxes 

Here we discuss generic properties of source, conditions of neutrino production and char- 
acteristics of neutrino fluxes. 



2.1 Properties of relativistic jets 

Hidden neutrino sources are associated to the core collapses of stars with masses < 
28M0. Models of these sources are based on extrapolation of the observed properties of 
GRB and models of observed jets. Recall that stars with mass > 28Mq and a fast- 
rotating core are widely believed to be the progenitors of the long-duration GRBs |20|. 



Evidences of highly relativistic jets, with bulk Lorentz factor Ff, ~ 10^ - 10^, have been 
found in recent GRB data [21|. 



A much larger number of core-collapses with masses of projenitors < 28M0 is 
believed to produce mildly relativistic, ~ lO"'^ - 10^, slow jets which do not break 



through the stellar envelope unlike the GRB jets [22|. The general picture is that materials 
from the central engine (a black hole or a highly magnetized neutron star created from 
the core-collapse) are emitted in lumps or shells with mildly relativistic speeds along the 
rotation axis of the star, thus forming a slow jet. Some initial shell or shells push out stellar 
material (with sub-relativistic or mildly-relativistic speed), thus make a cavity. Subsequent 
shells, with relativistic speed in the cavity, collide with each other due to variable outflow 
of accreting materials similar to the GRB internal shocks which take place well outside the 
stellar envelope (for reviews see e.g. [p^]). Each of these binary collisions produce shock 
waves (forward and reverse) in the colliding shells. The shell(s), initially ejected, can also 
produce a shock in the envelope (forward shock), and a reflected shock (reverse shock) in 
the shell(s). The shock in the envelope dissipates very quickly because of a higher density, 
but could also accelerate particles there. 



The subsequent binary colhsions between the shells take place mainly at the edge of 
the jet or at the inner border of the envelope. The shock waves that are generated in the 
leading and trailing shells, in a binary collision, are mildly relativistic, with the Lorentz 
factor 1 in the frame of the shells which, in turn, are moving at the Lorentz factor ~ 3 in 
the observer frame As these shock waves traverse the shells, magnetic field is generated 
from turbulence in the upstream and downstream regions of the shock front. The magnetic 
fields in these two regions are similar. 

The shocks in the hidden or hurried jets are optically thick to 7-rays as they are 
produced under the stellar envelope. High energy neutrino fluxes, however, are formed 
inside the slow jet due to interactions of shock-accelerated protons with matter, {pp— r 
collisions) and EM radiation (^7— collisions) of jets |8|. | , 

Following Refs. P, 0] we adopt a slow jet model with the following characteristics: 

• the total kinetic energy of the jet released over its duration is Ej ~ lO^^'^E'si.s ergs 

which is much smaller than the typical GRB jet energy; ^ — ^ 

• the jet duration is tj ~ lOtj^i s, which is the typical duration of the central engine's (jQ 
activity observed in GRBs; 

• the half angle of the jet, 9j ~ l/r^, is rather large and implies an isotropic-equivalent [\3 
total jet energy Ej^i^o » 2r^£'j; 

• with a variability time scale ~ 0.1i^^_i s and ~ lO'^'^Fb^o-S) the internal I ^ 
shocks (collision between two shells) take place at a radius rj « 2r^ct^ ~ 6.3 • C ^ 
10^°r2 g gt^^.i cm; 

• the width of the shocked shells is Arj ~ fj/'^l ~ 6.3 • 10^t„,_i cm and there can be O 
w tj/2ty ~ 50 consecutive collisions between shells which form shocks in the shells (jQ 
during the jet lifetime. | \ 

This picture is somewhat idealized since GRB data show more complex time structure 
than emission from identical shocked shells. Nevertheless, it captures the basic scenario of 
internal shocks. 

The pre-shock number density of particles in the jet calculated in the jet comoving 
frame (we denote the corresponding characteristics by "/") is 

/ / ^i,iso on 1 rv20 —3 -^51.5 -i \ 

n„ ~ n„ - — ~ 3.2 • 10 cm " —. ^ — . (2.1) 

(Here and below h = c = 1.) The post-shock number density is ~ An'p and the strength of 
magnetic field that forms due to turbulence in the shock region is 



B' « J327TeBn'j,mp ~ 6.3 • 10^ G J , (2.2) 



where es ~ 10 "^63,-2 is the fraction of shock energy that goes into creating magnetic 
field. The magnetic field drops down to the surrounding value of few Gauss (which is the 
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field from the "central engine", a magnetar or black hole) between successive shocks. The 
temperature of the thermal photons, created by shocked electrons which carry a fraction 
eg ~ 10~^ee,-i of the total energy, in the jet equals 

kT' « i - 4.3 keV [ 'r^^''^ \ • (2.3) 
2.2 Proton acceleration 

We consider a Fermi acceleration mechanism for protons in shocks. An acceleration time 
scale is proportional to the Larmor time scale and equals 

where Bcr = rn^/q ^ 1.488 x 10^° G is the critical magnetic field and (p ~ lOf/^i is the 
number of gyro-radii required to increase the particle energy by e-fold. We assume an 
acceleration spectrum N{Ep) oc E^"^ . Then the differential flux of protons, if they could 
escape freely from the jet, at a luminosity distance would be 

The maximum proton energy is determined by the shortest time scale for energy losses 
which is the synchrotron cooling time scale 

^9 1 /BerV 1 . ( KV' ( B' 



P 



In the jet frame this gives 



/Q ml Bcr -1/2 f B' , , 

and consequently, in the observer's frame -Ep^max = ^hE'^^^^ ^ 6.3 • 10"^ TeV. In turn, this 
energy determines the maximal energy of the neutrino spectrum. 

2.3 Neutrino fluxes from meson decays 

The rate of pp— interaction by shock-accelerated protons in the jet is given by 



~ 3 • 10^ s"^ 



( ^\ + 0.0073 ln2 



1 + 0.0548 In —4- +0.0073 In 



— — +u.uuYain" — — 

VTev; VTevy 



where cjpp is the total inelastic cross-section and the parameterization is valid for E'^ > 
10 GeV. The scattering rate of p^— interactions with thermal photons in the jet above the 



threshold of vr— production is given by 



JO e Jeth 

_ poo -/2 

~ ^ / , 10^ s-i. 

^ A;,mp/2i?;, exp[e7A;T'J - 1 

Here e'^ = — Pp cos 9) E'^/ nip is the photon energy evaluated in the proton's rest frame 
for the angle 6 between the directions of the proton and target photon. We used for 
simiplicity a constant, cross-section crp^(e^) ~ ctq ~ 200/x6, above a threshold photon 
energy e'^ ^ ^ 0.2 GeV for pion production in the rest frame of the proton. The 
threshold energy of protons which produce pions in the p'j— interactions equals according 
to O 

K,tk ^ ^ - 42 TeV. 

Thus, th.6 pp — scattGring rate doniinatcs below Ep^^ii ^ £'^^^r^ ~ 133 TeV, that is the 
whole energy range of interest. 

The fluxes of vr— and K— mesons from the pp— or pj— interaction at production can 
be calculated as 

dEp^p{Ep)Kpp/p^{Ep)t'^y^-^ — , (2.6) 

where Y^i^j^-j = Ep{dn^(^x^/dE„(^x-^) is the pion (kaon) yield function from the pp— or 
P7— interactions, and t'^^^ = Tj/Tf, is the dynamic or light crossing time. The quantity 
Kppip-y{E'p)t'^y^ is equivalent to the optical depth of the respective interactions and is very 
large for both the pp— and ^7— processes. For roughly constant values of Kpp/p^{E'p)t'^y^ 
the integral in (|2.6| ) corresponds to the fraction of the proton beam energy carried by vr"^ 
(/^+ ~ 17%) and {fK+ ~ 2%) mesons and follow the primary proton spectrum p^ . 
Thus the corresponding spectra at production can be estimated as 

^0 (F ■] f-K+{K+)Ej,\so 

^..,r...^i,^^K)} A^dlEl^j^^t^ HEp^^^jnmpY ^'-'^ 



following Eq. (^.4|). For the vr and K fluxes the fraction of the proton beam energy 
carried by the mesons are f^^- ~ 13% and fx- ~ 1%, respectively. Charm production and 



semi-leptonic decay contribute to neutrino flux at very high energies |25], which we ignore 
for the present discussion. 

We use the decay constants listed in Ref. to calculate the neutrino fluxes from the 
direct decay channels tt/K — t- , the chain decay channels vr /K — t- — t- i/^ z^g and the 

—7- TT^e^h'e decay channels by taking into account their energy losses (see Appendix A). 
The fluxes are plotted in Fig. |l]. Note that very high-energy neutrinos interact inelastically 
with nucleons in the stellar envelope and are subject to absorption. We discuss this issue 
in Sec. 2.7. 



The ratios of the Ve to and to fluxes can be fitted reasonably weU as 
<J)0 1 q 

'^^^ " if " (i^./GeV)i-65 [Q(^- - 10 + - ^^^^ 

+0.0265 [e(E^ - 45 GeV) + 9(10^ GeV - E^)] ; 
$Q 7 85 

+0.0265 [e(-E^ - 45 GeV) + 6(10^ GeV - E^)] , (2.8) 

where Q{x > 0) = 1 and 0(x < 0) = 0. The ratios decrease from the initial value because 
of a decreasing contribution to the i/g— flux from the chain vr — decays, and becomes 
constant at higher energy because of an approximately constant ratio between the fg— flux 
from — )• and f^— flux from — ;> /^^i^/i channels. 

The flux of muon neutrinos in the production region (the Vj ~ 6.3 • 10^^ cm) and for a 
jetted source at 10 Mpc can be parameterized, with a piece-wise power-law function, as 



GeV cm ^ s 



2 „-i 



' 1.64 ; 1 < E/GeV < 6, 
< 1.64 {E/6 GeV)-0-55 ; 6 < E/GeY < 600, 
, 0.13 (^/600 GeV)-i ; 600 < E/GeY < 3 x 10^ . 

This fit reproduces exactly the flux plotted in Fig. |l] below and above 6 GeV and 
600 GeV, respectively, and deviates by at most 30% from the numerical flux within 6- 
600 GeV. We use the numerical flux values for all our calculaions. 

The ratio of neutrino and antineutrino fluxes can be as large as ^ly/^p ~ 1.5 because of 
the preferential production of vr^ and over vr" and K~ . Indeed, in the pp-collisions the 
leading pion (which gives main contribution to neutrino flux) is positive. In pj— collisions 
there is a leading vr"*", although the secondary neutron with roughly 85% of the initial 
proton's energy will interact further and produce a vr^. 

This effect is present in the atmospheric neutrino fluxes, however the ratio of neutrino 
and antineutrino fluxes and details of the energy dependence is different in our situation. 
The original flavor content is thus, 

e{E) : 1 : and e{E) : 1 : , (2.9) 

where e{E) and €{E) decrease from 0.5 at low energies down to 10^^ at high energies. 

The charged vr, K mesons and muons can also be accelerated in the shocks |2^ similarly 
to the protons, with a maximum energy E.,^-K,ma.x = {'mTT-K / mp)'^ Ep^y^ax following Eq. (^). 
Although a detail study is lacking, it might be possible that vr, K and fj, decay without 
suffering severe energy losses. In such a scenario, the flux ratio 1:2:0 will be maintained 
to high energies. In reality one can expect some intermediate situation between this flavor 



content and the one in Eq. (2.9) without acceleration of mesons. In what follows we will 
present computations for both extreme flavor contents. 

The parent mesons and muons are isotropically disributed in the shocked region be- 
cause of their small Larmor radius in the large magnetic fleld. Therefore the neutrinos are 
also isotropically produced in the comoving frame. However in the observer frame they are 



emitted mostly along the jet direction similar to the relativistic beaming effect for photons. 
Then neutrinos, produced in the low-density jet, propagate through the stellar envelope. 




Figure 1: Neutrino and antineutrino fluxes (thick curves) at production from a hidden source 
at 10 Mpc. The red (purple) thick solid curve corresponds to {ve) flux while the red (purple) 
thick dashed curve corresponds to (Pg) flux. The labeled thin dot-dashed curves correspond to 
different components of the i/f^ flux. 



2.4 Stellar envelope: density profile 

The pre-explosion stars of the GRBs and the associated type Ic supernovae are widely 
believed to be He stars with radius i?^ ~ 10^^ cm. Supernovae of type 11 and lb are 
thought to be explosions of larger stars with radius i?^ w 3 x 10^^ cm, and have a H 
envelope on the He core. The envelopes usually have density profile p(r) oc r~^, with 
A; = 1.5 - 3 |2^. Following Ref. ||l^ we adopt two stellar profiles, model [A] and [B] 
parameterized as 

[A] p{r) = 3.3 • 10~6 g cm-3 - l 

\ r 

[B] p(r) = 2.8 • 10-5 g cm-3 | ^^^1'^'''': ^0 < \< ^ (2.10) 

The radius of inner border of the envelope equals Rq ~ Vj ~ 6.3 • 10^*^ cm. According to 
Eq. ( 2.10| ), the density at Rq is p^ 0.33 g cm^'^ but can vary between (0.1-1) g cm^^. 
The corresponding density of electrons is ng = pqNaYs ~ 10^'^ cm~^ for Yf, = 0.5. Thus, 
neutrinos are produced in a low density region of jet and propagate through a higher density 
envelope. We assume that the boundary between the jet front and the envelope is sharp 
(much smaller than any oscillation scale). Neutrinos are produced in a low density region 
of jet. Typical size of this region is about ~ 6.3 • 10^ cm, that is, about 1/10 of a jet 
size. 

Notice that in general, all characteristics of a star change in time over jet duration. 
In particular, the jet parameters such as Tf, and Ar^ are subject to variation over time. 



The jet radius rj is expected to increase slowly with time as the jet burrows through the 
envelope and eventually chokes or successfully breaks out, the parameters of the envelope 
Rq rj and no change with time. 



2.5 Neutrino fluxes and densities 

According to Fig. |l| the flux integrated over the energy above 1 GeV, which gives approxi- 
mately the total flux of neutrinos at rj, is -F(> IGeV) w 10'^'^ cm~^ s~^. Then the number 
density of neutrinos in the source equals ~ 10^^ cm~^. Notice that this is a densitity in 
the reference frame of the jet where neutrino emission is approximately isotropic. Therefore 
the potential due to — scattering equals /i = yplGpr)!^. It should be compared with 
the vacuum frequency uj = ^m^^j^E' . We find that ~ a; at the energy E' ~ 10^ GeV, 
which corresponds to ii^ ~ 3 • 10^ GeV in the observer frame. Consequently, for energies 
of interest the neutrino-neutrino effective potential is too small to induce the collective 
neutrino effects. Some collective effect may show up at > 1 PeV. 

2.6 Coherence and averaging 

Let us esimate the size of the wave packets of the produced neutrinos |18]. Pions mainly 



decay near shock fronts in the regions with strong magnetic field and increased density. 
Pions undergo collisions with surrounding photons and gas. The latter (gas) dominates 
at low energies (£^ < 4 • 10^) GeV. The mean free path with respect to the collisions 
can be estimated as X^oi = (10^, 3 • 10^, 5 • 10^) cm for energies E' = (0.1, 1, 10) TeV, 
correspondingly (in the rest frame of jet bulk). For neutrino emission in the forward 
direction the size of the wave packet is then |jl^ ax = ^coi/Itd where 7,^ is the Lorentz 
factor of pion. This gives ax = (10^, 3, 5 • 10~^) cm for energies E = (0.1, 1, 10) TeV, 
correspondingly. However, very strong kinematical shortening of the wave packets occur 



in presence of the magnetic field due to bending of pion trajectory. According to |18], the 
size of the wave packet is given by 



.. = 3.5.10-- cm {-^\ . (2.11) 



For Pf, = 3, S = 10^ G and neutrino energies E = (10^, 10^, 10^) GeV we obtain from 
this formula ax = (10"^^, 3.5 • 10"^^, 10"^^) cm correspondingly. 

Separation of the wave packets of different mass states on the way L is given by 
ds = LA.m? / {2E'^). For typical distance 10 Mpc and Am|;^ we find 

/I TeV\^ 
ds = 0.045 cm ( — ^ J . 

Thus, ds ^ ax, and therefore the wave packets are well separated at the detection site. 
This means that the coherence is lost in configuration space in the course of neutrino 
propagation. 

The coherence will not be restored at the detection. Indeed, the separated packets 
interact coherenly, if the detector has long memory (time interval of coherent detection). 



In this case the packets interactions would interfere and produce oscillatory pattern in the 
energy scale. The period of this oscillatory pattern, that is, the energy interval over which 
the phase changes by 27r is given by /S.Et = El,^/L, where li, is the oscillation length. For 
L = 10 Mpc and Am^]^ we obtain 

Ai?T = 3 • 10-3 eV (^p^) . (2.12) 

In practice it is not possible to determine the neutrino energy with such an accuracy. 
According to Eq. (|2l^) at 1 TeV one needs to have the energy resolution AE/E < 10 to 



see non-averaged oscillation effect, whereas one may achieve AE/E ~ 0.1. Thus, incoherent 
fluxes of mass eigenstates arrive at the surface of the Earth and interact in a detector. 

Let us estimate the coherence length (the distance over which the wave packets are 
completely separated): Lcoh = crx2E'^ / Am? . Using ax from ( 2.11| ) we obtain for Am'^i = 
2.4 • 10-3 eV^ 

\3 B J \l TeV 

Notice that Lcoh oc \fE and for E > 200 GeV, Lcoh becomes larger than the size of star. 

Some small part of pions decays in jet between shocks where the magnetic field is much 
smaller. For 5 = 10 and E = \ TeV we obtain Ox = 3.5-10-^° cm and Lcoh = 1.4-10^''' cm. 
Similar situation is for the 2-body decay of K-mesons. 

For muon decays (and also for 3-body K— decays) in the magnetic field, size of the 
neutrino wave packets is much larger p8[ ]: 

/lO^ G\ /I TeV^ ^ 



ax = 1.7- 10-^^ cm 
For Arrig^ the coherence length equals 



3 / V B \ E 



Lcoh = 7 -10^^ crni^j f — ^ 

and it does not depend on E. Although this length is about 2-3 orders larger than for 
pions, conclusions about loss of coherence are the same. 



2.7 Effect of inelastic interactions 



Evolution of the flavor neutrino states uj = (fg, i^n, Vr) is described by the equation: 



i^ = {H + Hint)uf, (2.13) 

where H is the standard Hermitian part which includes the vacuum and refraction terms 
(the real part of scattering amplitudes, A^'^) and Hint describes inelastic interactions (the 
imaginary parts of scattering amplitudes, A^"^). Since the interactions are flavor diagonal 
in the lowest order of perturbation theory, we have 



Hint = --;:diag{Tc,T^,Tr). 



For E > 3 ■ 10^ GeV (which corresponds to 0.5r = A^"^ ~ A^'^ for scattering on nu- 
cleons) the inelastic interaction can not be neglected. However, for such high energies 
in the first approximation the inelastic amplitudes for all neutrino species are the same: 
Fg ~ ~ Ft- = r. Indeed, the difference of masses of the charged leptons as well as 
the inelastic scattering on electrons can be neglected. Consequently, the inelastic part of 
the Hamiltonian becomes proportional to the unit matrix: Hint = Ii^ this case the 

inelastic interactions and oscillations factor out. The former does not influence the oscil- 
lation pattern. Indeed, we can define new flavor wave functions as Uf = e^^/'^uj. Inserting 
this relation into ( 2.13| ) we obtain 



i—r- = Hvt 

dt ^ 

without the inelastic part. The factor due to inelastic interactions then appears in the 
probabilities: P = e~^*P. This factor describes absorption of neutrinos or scattering which 
lead to departure of neutrinos from the coherent state. The effect of neutrino absorption in 
the charged current processes can be included in the energy spectrum of neutrinos. Decays 
of n and r produced in the CC interactions of neutrinos will generate secondary neutrinos. 
Due to falling down spectrum with energy, the contribution of these neutrinos to the total 
flux is small. In the case of neutral current neutrino interactions one should consider 
oscillations of the scattered neutrinos. In what follows we will neglect this scattering since 
the event rate of neutrinos with E > ?> ■ 10'^ GeV is small even for the nearest plausible 
source. 



3. Conversion inside the star: general consideration 
3.1 Propagation basis 

For the mixing matrix in vacuum, Vf = C/pMNS^^; we use the standard parameterization 

UPMNS = U23{923)rsUi3{9i3)Ui2iei2), (3.1) 

where Uij = U{6ij) is the matrix of rotation in the zj— plane over the angle 9ij and F^ = 
diag{0,0,e^^). It is convenient to consider dynamics of the conversion in the so-called 
propagation basis, v'j, defined as 

l^f = U23T5l^'f , l^'f = {ve,v'^,v'^) ■ (3.2) 

In this basis the CP-phase is eliminated, so that the dynamics of conversion is not affected 
by 5. The CP-violation appears only in projection of the initial and final states onto the 
propagation basis. We have explicitly 

10 

U23TS =10 C23 S23e'^ I , (3.3) 



>0 -S23 C23e 



iS 



where S23 = sin ^23 , C23 = cos ^23 ■ According to ( ^.iD and (|3.2[) the vacuum mixing matrix 
in the propagation basis is given by v'r = U'v, where 



U' = UisUu 



C13C12 C13S12 Sl3 
-S12 C12 

-S13C12 -•S13S12 Ci3, 



(3.4) 



Here C12 = cos0i2, S12 = sin0i2, etc. It is straightforward to show that the mixing matrix 
in matter in the propagation basis, U^, up to additional smah 2-3 rotation has the same 
form as in (3.4) with mixing angles in matter: 612 



am 
'^12' 



'13 



'13- 



„m „m 
'-13''12 



•'12 
,r7i m 
'13''12 



'-13*12 



"-12 
-,m m 
'13*12 



'13 



'13 ^ 



(3.5) 



The level crossing scheme (see Fig. |2|) is similar to the one for low energy supernova 
neutrinos (see, e.g., ||3^ ). 



E,„ (eV) 




n,, (cm ) 

jq23 



[t] 



Figure 2: The level crossing scheme for neutrinos in stellar envelope in the case of normal mass 
hierarchy. The solid lines correspond to the eigenvalues of the effective Hamiltonian as functions 
of the electron number density. The dashed lines correspond to the energy levels of the flavor 
states i^e, v^i and Vr' in the propagation basis. The vertical dot-dashed lines indicate the L- and 
H- resonance densities. Note that the positions of the resonances are inversely proportional to the 
neutrino energy E (plotted here for £ = 2 TeV). We used sin^ 2^13 = 0.08. 
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3.2 General expressions for probabilities 

Decays of vr— , /i— and K— mesons produce incoherent neutrino and antineutrino fluxes of 
definite flavors: v^, f/^, etc.. Since the mass states are the eigenstates of propagation in 
vacuum and coherence between them is lost due to separation of the corresponding wave 
packets (see Sec. 2.6), we can write for the Ua — ^ I'p conversion probability at the Earth 

^'(^^a ^ i^/?) = J^P*(i^a ^ i^i)|C//3iP, (3.6) 

i 



- 13 - 



where — I'i) is the flavor-to-mass conversion probabihty inside the star. This ex- 



pression differs from the one used in Ref. |16|, and detailed explanation of the difference is 
given in the Appendix B. If neutrinos cross the Earth the formula (|3.6| ) is modified as 



(3.7) 



where PE{^i — ^ ^p) is the probability of {vi — >• vp) oscillations in the matter of the Earth. 
The probability Pe can substantially deviate from |C/^ip at energies E < 10 GeV. For 
higher energies matter suppresses oscillations inside the Earth, so that the probability is 
reduced to the one in Eq. (^.6|). We give explicit expressions for the probabilities in the 
Earth for constant density in Appendix C. 

Inside a star the flavor conversion occurs first in jet and then in envelope. Therefore one 
should take into account oscillations in the production region of jet and perform integration 
over this region. Then the total conversion probability inside the star equals 

2 



r 



dx 



V JO 



where Aj(.t{va — ^ ^5)(^) is the amplitude of probability of Va transition to the flavor state 
inside jet: between the production point x and the inner border of the envelope. The 
integration is performed over the production region and in what follows we will use the 
one-dimensional (ID) integration for simplicity. Note that ID consideration is valid if we 
consider neutrinos from a given source, although in this case the density proflle can differ 
from the one in radial direction (recall that jet has rather large cone angle). If Ot^ 



> 



'jet 



IS 



the angle between the radial (jet) direction and direction to observer, the radial distance 
equals r = x cos 6tr, where x is the distance along the trajectory. Then the density profile 
which neutrino experiences is n{xcos9tr) = nQ{x cos Otr / f'o)~^ ■ In this case the density 
profile is flatter and the evolution will be more adiabatic. For neutrinos produced by all 
hidden sources (the diffuse flux) we need to perform also integration over cos Otr- 
In the propagation basis the probability can be written as 

2 



r,. 



dx 



X]X^(^23r5)a/3 • Ajet{v'is v's^){x) • Aenv{v'^ Vi) 

i /3 



(3.^ 



Here Ajeti^'j^ v'^ is the amplitude of transition inside jet in the propagation basis. The 
probability averaged over the production region equals 



P*{Va Vi) 



1 

Jo 



dx Sin 



where Sia = A^:{i'a Vi). The matrix of total flavor-to-flavor probabilities equals 



P 



tot 



P*Pe, 



(3.9) 
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or if there is no Earth matter effect: 



tot 



p.-p> 



PMNSi 



(3.10) 
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where PpMNS = II \{UpMNs)ij\'^\\ is the matrix of moduH square of the mixing matrix 
elements in vacuum. 

According to ( |3.8D the amphtudes of transitions inside the star is given by 

In matrix form we have 

5"* = {U2Z^s)SjetSenv (3-11) 

Here 

{Sjet)il3 = Ajetil^'is f^) and {Senv)i( = ^eni>(l^5 l^i) 

are the evolution matrices inside jet and inside an envelope in the propagation basis cor- 
respondingly. 

For E > 2 ■ 10^ GeV or / and smaller production region the phase of oscillations inside 
jet is small and these oscillations can be neglected. Furthermore, as we will show in Sec. 4.5, 
the effect of oscillations inside jet is zero or small and so in the first approximation can be 
neglected. Therefore we will first consider oscillations inside the star ignoring oscillations 
inside jet. In this case Sjet = I and S** « {U23Ts)Senv 

3.3 Matter affected range 

In Fig. ^ we show the conversion probabilities in different channels for two different initial 
densities in the envelope. There are two key energies in the problem: the two resonance 
energies which correspond to maximal electron density in the envelope no (density at the 
border between jet and envelope) ^: 

where Vq = ^/^GpnQ. For uq = 10^^ cm^^ we find 

Ej^ = 1.3 GeV, E^ = 75 GeV. 

These energies determine the borders of energy regions with different dynamics of flavor 
conversion (see Fig. ^) inside the star: 

1. E < E^- the vacuum oscillations (VO) region: VO dominate for both 1-2 and 1-3 
modes. 

2. E^ ^ E < E^ - the intermediate energy range, (7 — 70) GeV. The corresponding 
neutrinos are produced between the two resonances. The 1-3 mixing and split lead 
to averaged vacuum oscillations (the mass eigenstate fs decouples from the rest of 
system), whereas for 1-2 mixing and mass split the matter effects dominate. 

3. E^ ^ E < Ena, where Ena is the energy of strong adiabaticity breaking. Neutrinos 
are produced above the 1-3 resonance in the density scale; here matter effects are 
important for both mixings. For the power dependence of the spectrum we have 
Ena ~ i?*Amii/47r - W^E^. 



""^This is 21/— definition, in fact there is some small shift of the 1-2 resonance due to 1-3 mixing. 



4. E > Ena - matter suppresses oscillations inside the star, here the flavor conversion 
is due to oscillations in vacuum from the surface of a star to the earth. Loss of 
coherence leads to the averaged oscillation result. 




0.6 




"10° 10' 102 10-^ 10"* 10^ 

E [GeV] 

Figure 3: Probabilities of the v^. — vp (top panel) and vp (bottom panel) transitions for 

the density profile A with two different inner densities no = 10^"^ cm~^ (red lines) and no = 2 • lO^'^ 
cm~'^ (blue lines). We used sin^ 20i3 = 0.08, sin^ (?23 ~ 0.5 5cp — and normal mass hierarchy. 

The matter affected range spans over 5 orders of magnitude from about 0.3 GeV to 
3 • 10^ GeV. With increase of no the range expands in both directions. ^ According to 
Fig. ^ generic features of the energy dependences of the probabilities can be summarized 
as follows. 

• Plateau in the intermediate range above the vacuum-to-matter transition region at 
(0.3 - 3) GeV. 

• Peak or dip (depending on channel) at > E^ due to the 1-3 resonance. 

• Wiggles at £' > E^ in the nonadiabatic edge, above the peak (dip). As we will show 
the wiggles are due to interference induced by adiabaticity breaking. 

^Here we discuss whole energy affected range and discuss relevant physics in spite of the fact that the 
present detectors have rather high energy threshold E > 10"^ GeV. Future experimental developments may, 
in principle, substantially reduce this threshold. 



• The wiggles are different in channels with initial Ve and u^. According to Fig. ^ the 
wiggles are small in Vf, — )• i/^ channels (upper panel), and large in — ?■ 1^/3 channels. 

• The wiggles modulate the non-adiabatic edge of the energy profile of the effect. 

Note that the Ue probability as function of energy can be viewed as two dips 

in the energy profile, one embedded into another, which corresponds to crossing of two 
resonances: due to the 1-2 and 1-3 mixings. If the 1-3 mixing is small or zero, the 1-3 dip 
is absent and the plateau extends to higher energies. With further increase of energies the 
probabilities will converge without wiggles to the asymptotic values as in the case of solar 
neutrinos. In what follows we will explain these features. 

3.4 Asymptotic values of probabilities 

Outside the matter affected region the probabilities converge to the asymptotic values 
which are given by the averaged vacuum oscillation probabilities. Indeed, for E <^ 
the matter effects can be neglected and we deal with vacuum oscillations inside the star, 
(essentially the flavor state is projected onto mass eigenstates in the production region and 
then these eigenstates propagate inside the star without changes), so 

P.{Ua ^ yi) = \Uai?. (3.12) 

and therefore Eq. (|3.6| ) becomes 

i 

The same result ( plD , ( |3l^ ) are valid for very high energies: E ^ E^. In this region 
matter suppresses oscillations inside the star. In other terms, for very high energies the 
adiabaticity is strongly broken and no flavor conversion occurs, P.^ = I. Here we project 
onto mass eigenstates at the surface of the star, which is equivalent to the coherence 
lost. Again P*(t'a — )• Ui) equals ( 3.12| ) and the total probability is given by (|3.13 ). The 



probabilities reach the asymptotic at ~ 10^ GeV and converge to the same asymptotic 
values below 1 GeV (in this sence the picture is symmetric) as in Fig. |3|. 

Since for antineutrinos in vacuum Uai C/*j, the result (|3.13|) is also valid for antineu- 
trinos: P{i>a T^p) = P{^a ^p)- There is no CP-asymmetry, although the probabilities 
depend on the phase 5. Since the asymptotics does not depend on Am^, for the inverted 
mass hierarchy we have again the same result ( p^ . 



3.5 Adiabatic conversion 

At low energies {E < E^ ~ 100 GeV) the adiabaticity condition is satisfied and we can 
use results in the adiabatic approximation. Let i'^ be the eigenstates of the Hamiltonian 
in the initial moment (at the inner border of the envelope). Then the flavor states can 
be represented in terms of as i^f = U^v^. The adiabatic evolution means that the 
transitions between the eigenstates can be neglected and therefore — )• v, or explicitly 

(3.14) 




Consequently, inside the star the flavor state Ua evolves as Va ^ai^i^ ^^'^ therefore the 
adiabatic probability is given by 



rO |2 
ail • 



Then the flavor transition probability at the Earth equals 

Let us describe the adiabatic transitions in the propagation basis. Now Uj = U^i'^, 
where is the mixing matrix in matter in the intial moment (at the density no). Ac- 
cording to Eq. ( |3.2|) , the flavor states can be written as Uf = {U2z^5)SjetU'^v^. Therefore 



according to ( 3.14 ) Svf = {112^^8) SjetU'^v-, and 

= {u\S\uf) = {U2zTs)S,etUt 
Neglecting oscillations inside jet we have Sjet = I and 

= {v\S\vf) = {U2zTs)Ut (3.15) 



The expression in ( 3.15| ) coincides with the mixing matrix in matter in the initial moment: 

The adiabaticity condition can be satisfied at low energies, in particular, in the inter- 
mediate region. This region corresponds to neutrino production above the 1-2 resonance 
density. Far above the 1-2 resonance 6^^ ~ vr/2. Therefore using Eq. ( ^.5]) we obtain the 
mixing matrix in the initial state 
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-1 











*13 


Jd 

^13 



U'^^\-l 1, (3.16) 
where 0^3 and are the mixing parameters in the initial moment. Inserting this matrix 



into Eq. ( 3.15 ) we obtain the matrix of amplitudes 




^0 „o 

S* -Ki=\ -C23 -s\^S2ze'^ c\^S2^e'^ | . (3.17) 

-S?3C23e*'^ -C?3C23e*^ 

Consequently, the matrix of probabilities inside the star equals 

U <^13 *13 \ 

C23 ^13 ^23 "^13 '^23 I ' (3.18) 
„2 2 2 2 2 / 
*23 *13 '^23 '^13 '^23 / 

where the rows correspond to the initial flavor states v^^v^,, Vr and the columns correspond 
to the flnal mass states. According to ( ^.18 ), the probabilities do not depend on 5 and this 



is the consequence of the adiabaic propagation at the energies above the 1-2 resonance. 
So, any dependence on 5 is manifestation of the adiabaticity violation or/and closeness to 



the 1-2 resonance. The flavor transition probabilities at the detector are then given by 
Eq. (3.10). Some dependence on 5 follows from dependence of the projections of the mass 
states to flavor states back at the detection, i.e., from factors IUpmnsI"^ as we will discuss 
in Sec. 4.6. 



Let us consider some particular cases of Eq. ( 3.1S| ). 
1. In the intermediate range (plateau) we have 6*5*3 ^ probability matrix 

becomes 



'^13 ^13 



-P* ~ C23 S13S23 C.\3^23 • (3.19) 



i2 ^2 ^2 ^2 ^2 
'23 *13''23 ''13'^23 • 



2. Above the H-resonance, E > E^, we have ^^3 ~ '7r/2, so that c^^ ~ and 55*3 ~ 1, and 
consequently, 





t/™ « I -1 



-1 

For large enough ^13 the adiabaticity is fulfilled, and using the adiabatic result ( t3.15|) we 
obtain 



which leads to 









-C23 


-S23e 


S23 


-C23e 
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1 "^23 


*23 




r2 


V ^23 


•-23 



i6 
iS 

P*^ I 4, s^3 oj . (3.20) 
Then, according to general formula ( 3.10| ) the flavor probability matrix equals 

/ \Ue3\' M \Ur3\' \ 

P« ci3|C/el|2 + si3|C/e2p 4,\U^i\^ + sl^\U^2\^ C^Urll^ + ■43\Ur2\' ] ■ (3-21) 
\sUUel\'+43\Ue2\' S^U.ll' + 43\U,2\' ^Is I ^^rl P + 4 1 [7.2 1 V 

3. For antineutrinos, the 1-2 mixing is suppressed at > Ej^, so that ~ 0. Therefore, 
according to (|3.5D we obtain 



U!2^ I 1 



-13 *13 



^^3 C^i3, 



The matrix of amplitudes in the flavor basis equals 5* = (C/23r_5)C/^, or explicitly 



C?3 .?3 

~ I -S?3S23e-*'^ C23 C?3S23e-'^ 

-S?3C23e-*'^ S23 C?3C23e-*^ 

Consequently, the matrix of probabilities inside the star is 

02 n 02 

•^13 '^13 

„0 2^2 „2 „0 2„2 



2 2 2 2 2 ^0 92^ 

*13 ''23 '-23 '-13 '^23 • [^O.ZZ) 



2 2 2 2 2 
. *13 '-23 *23 '-13 '-23 • 



4. For the inverted mass hierarchy (IH) above the 1-2 resonance we have the same initial 
mixing matrix as for the normal mass hierarchy (NH) given in Eq. ( p^ ), the same b- 
matrix (|3.17| ) and the same matrix of probabihties ( 3.1^ ). 



In the intermediate energy range (below 1-3 resonance) the 1-3 mixing equals approx- 
imately the vacuum mixing as in NH case and therefore the probability matrix is as in 
Eq. ( pQl ). The difference appears at higher energies, since there the matter effect on 1-3 



mixing is different for the normal and inverted hierarchies. In particular, above the 1-3 
resonance we have in the case of IH: 9^^ « 0, c?^ ~ 1, s?, ~ and 
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^23 / 





P* « C23 si, . (3.23) 



Then the matrix of flavor transition probabilities in case of IH becomes 

/ M |f^M2|' If^r2p 

P« ci3|[/elP + 4|C/e3p 4 1 C/^l | ^ + S^g | [7^3 P cl,\Url\^ + S^Ursl 
\si3|C/elP + ci3|C/e3p 1 C/^i | ^ + 0^3 1 C/^3 P | ^ + 0^3 1 1 

as compared to the NH case ( |3.21| ). 

For antineutrinos, the situation is similar: in the intermediate region the results coin- 
cide with those for NH, whereas above the H-resonance energy the results change. Now this 
is the resonance channel and c^, ~ 0, ~ 1. Consequently, the matrix of probabilities of 
the flavor to mass transitions becomes: 



(3.24) 



3.6 Adiabaticity breaking 

At high energies, not far above the H-resonance for the profiles A and B, the adiabaticity 













*23 
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r2 
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•-23 


*23 



is broken and in contrast to Eq. (3.14) the transitions between the eigenstates of propaga- 
tion in matter occur. The adiabaticity is broken mainly in resonance regions, so that in 
resonances jumps from one eigenstate to another occur. 

We will consider evolution in the so-called factorization approximation when two level 
crossings (1-3 and 1-2) occur independently one after another and do not influence each 
other (see Fig. |2|). Probability of crossing of two resonances equals the product probabilities 
of crossings of individual resonances. 

The "jump" transition probability in each crossing is well described by the double 



exponent formula |33] 



g27rA„(.ij(l— sin^ 9) 

Pjump = ^2ttX„uj _ I ' (3.25) 

where uj = Am'^/2E and A„, = n/{dn/dr) is the scale of the electron density change. For 1- 
3 level crossing the jump probability Ph = Pjumpi^iT^'ii, Ois) and for the 1-2 level crossing 
the jump probability ^^(^"^21' ^12) (see Fig. |2|). In Fig. Q we show the jump probabilities 
computed for different neutrino parameters. 



Let us consider Vf. production with energies substantially larger than the resonance 
energy, so that in initial state Vf. ~ vzm- In factorization approximation we have 



P,{ve^vi) = PHPL, PM^U2) = Ph{1-Pl), P.{ue^^^) = l-PH- (3.26) 
Then the {v^ — >• fe) probability equals 

P{Ve ^ Ve) = (1 - P//)|t/e3|' + Ph{1 " PL)|C/e2|' + P^Pl | f/el |'. 

It can be rewritten as 

P{Ve ^ ye) = \Uei? + PH{\Ue2? " We^?) + PHPL{\Uel? " \Ue2?)- 

For other channels we have 

P{ve ^ v^.) = (1 - PhW^z? + Ph{1 - Pl)\U^,2? + PhPl\U^i\\ (3.27) 

and similarly for transition to i^r- 

The jump probabilities have the following asymptotics for £' — t- oo: 

PH^l-\Ue3\\ Pl^y^^^. (3.28) 

Indeed, the limit -E — )• oo corresponds to absense of flavor transformation, i.e., — )• fe- 
If the initial state is fe = fsm, the result for Ph is obtained in the 2i^— approximation 
neglecting matter effect on the 1-3 mixing in the intermediate point between the two 
resonances. Then the expression for P/, can be found from the condition that the Su— 
probability P*(z^e — ^ i^i) = l^^eiP (since —?■ Ug), and therefore according to ( 3.26| ): 



Pf, = iC/eil^/Pff- Finally, using asymptotic expression for Ph we arrive at the result shown 



in Eq. (|^ 



It is interesting to note that the adiabaticity is broken quite similarly in the same 
energy range for both resonances. This feature is the consequence of the fact that for 
sin^ 2013 ~ 0.08: sin^ 26*13 Am|i ~ sin^ 2612 Aml^. 

Due to breaking of adiabaticity new interference effects emerge which lead to oscillatory 
behavior ("wiggles") of the conversion probabilities (see Fig. |^. There are two types of 
the interference effects which are related to the adiabaticity violation in the H-resonance 
and L-resonance. We will call the results of these interferences the H-wiggles and the L- 
wiggles correspondingly. In general both L- and H- wiggles are present in all the channels 
simultaneously, however, as we will show in a given channel one type of wiggles dominates. 

3.7 H-wiggles 

The H-wiggles are related to mixing of the eigenstates i'2m and 1/3^ in a given flavor state 
and to the adiabaticity violation in the H-resonance. Above the H-resonance substantial 
admixture of i^2m and 1^3^ appear in the fg— state, and therefore the H-wiggles dominate 
in channels. 

Let us neglect 1-2 mixing and consider transition — >• z^2m. in the region of H-resonance, 
i.e. compute the probability to find z^2m. at densities below the H-resonance but above 
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Figure 4: The jump probabilities in the H-resonance Ph (black lines) and in the L-resonance 
Pl (red line) as functions of neutrino energy. Ph is shown for three different values of sin^6'i3. 
We use the profile A and neutrino parameters: Am^^ — 8 • 10^^ eV^, Am§j = 2.4 • lO^'^ eV^, 
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the L-resonance. The problem is reduced to 2i/-problem according to the factorization 
approximation. In the H-resonance there are two channes of transition z^g — ^ '^2m- 
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where the first arrow denotes projection of Vf. state on the corresponding eigenstate. It is 
the interference of the corresponding amplitudes that leads to the H-wiggles. 
Let us quantify the interference picture. At the production point 



Ve ~ cos V2rn + sin Vj,m-, 



(3.29) 



where ^J*, is the mixing angle in matter at the production point: 



sin^ 26l?3 



sin^ le 



13 



cos2 2013(1 - EjE^ )2 + sin^ 2^13 ' 



Above the resonance, E > E^^ the angle ^13 ■nj'l and the admixture of V2ra^ given by 
cos 05^3, decreases. According to Eq. ( 3.2g| ), below the H-resonance the total amplitude of 
probability to find v^m is given by 



■?,^22 + sin Q%.A, 



cos Vl-^J\22 



13^32, 
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where A32 is the amphtude of v^m — ^ i^2m transition in the H-resonance, 1^32] = VPh, and 
1^22! = \/l — Ph is the amphtude of survival probabihty. Then the probabihty to find i^2m 
equals 



P{l^e ^ 1^2m) = COS2 0?3(1 - Ph) + sill^ O^Ph " sin /FV(W^ COS , (3.30) 

where 

(f)H = Arg{A22Al2). 

The last term in Eq. ( |3.30| ) is the interference term of two amplitudes which is responsible 
for the wiggles. The amplitude of wiggles is given by 



W = sm2e^^^y/PH{1-PH). 



Using the formula ( 3.25| ) for Ph and for sin^ ^13 = 0.02 we find the following values: 



W = (0.05, 0.035, 0.025, 0.013), for E = (200, 300, 500, 1000) GeV respectively; in 
agreement with results of the numerical computations in Fig. ^. 

Let us estimate the phase (pn- Note that transition between the eigenstates 1^3^ — ^ ^^2™ 
is driven by the derivative 9^ and occurs mainly in the resonance region it Arjj. Here 
'^R = fR{E) is determined from the resonance condition, \/2GFne{rR) = Am'^i/{2E), and 
the half-size of the resonance region is given by 

Tl 

ArR = A„ tan 26*13 = , , , tan 26'i3. 

an/ar 

For the power dependence of the density on distance, 

\ -k 

the scale factor equals A„ = r/k. Therefore Arji = r/j tan 2^13/^:. 

The whole evolution of the phase <j) can be devided into three parts: (i) from Rq 
to rpt — Ava, where the system evolves adiabatically and transitions ly^m ^ T^2m can be 
neglected; (ii) from — Ar^ to + Ar^, where the transitions 2/3^ ^ i^2m occur; and 
(iii) from vr + Ar^ to i?*, where the adiabaticity is restored and states usm, i^2m evolve 
independently again. For moderate adiabaticity breaking Ar^ < Ar/j; for high energies 
the adiabaticity is broken even beyond the resonance layer. 

The total amplitudes j422 and ^432 introduced above can be written as j422 = ^22^22^22 
and A32 = ^22^32^33 since ^432 7^ only in the region (ii). Note that the amplitudes in 
the region (ii) and consequently phases coincide for both yl22 and ^32. Correspondingly, 
the total phases (/>2 = arg{A22) and (ps = arg{A32) equal 

CP2 = ct>f + <t>? + <t^'^ , ct>, = <pf + 4^ + 0f ) , (3.31) 

where 0^*^ and 0^*"'' (j = 2, 3) are the adiabatic phases in the regions (i) and (iii) cor- 
respondingly, whereas 02*^ = ^^^5(^22^) ^^"^ "^s*^ = '^^5(^32^) are non-adiabatic phases. 



Note that the last terms in Eqs. ( |3.31| ) (contributions from the region (iii)) are identical. 
Thus, the phase difference we are looking for equals: 



(ii) A^^)*^ 



23 



where Acf)^^^ = (jr^ — (jj^ is the adiabatic phase difference acquired in the region (i): 

frR{E)-Ara{E) fr Am^T^ /Aw^X ^ 

A0W«y^ dxMv{x)-cos29,3^^ +sin22%(^^j . (3.32) 

For power-law density profile the upper limit of integration equals rji{E){l — tanlOi^/k). 

Few comments are in order. The phase (pn does not depend on evolution above (out- 
side) the resonance (region (iii)) and therefore it does not depend on R^,. The phase depends 
on rj via the adiabatic phase contribution. The non-adiabatic contribution Ar(7[^^^^*] < 
27r (or much smaller) since the adiabaticity violation corresponds to the width of the reso- 
nance layer to be smaller than the oscillation length. To a good approximation the phase in 
region (ii) can be taken into account by extending the integral in Eq. ( [3.32 ) to this region, 
that is, taking rji{E) + Ara{E) as the upper limit of integration. 

If y ^ Am? /{2E), the adiabatic phase can be estimated in the following way: 

A(^(*) « / V{x)dx - —=r cos 2^13 [vr - Avr - Rq] , (3.33) 

where we have taken Ar^ = ArR. This allows us also to estimate the period of wiggles in 
the energy scale. The first term in ( |3.33| ) only weakly depends on energy. This is because 
V{x) is a steeply decreasing function of x (the power /c ~ 3) and the integral is given mainly 
by the lower limit V{Ro), which is independent of E. The upper limit depends on E but 
its contribution to Acj)^'^^ is smaller than the second term in ( 3.33D and therefore will not 



change the period of wiggles substantially. Hence it is enough to study energy dependence 
of the second term which can be rewritten as 

(padiE) « — [rR{E) - AvRiE) - Rq] . (3.34) 

This dependence on energy agrees well with periods of wiggles obtained from exact numer- 
ical computations (Fig. ^). The key point that the phase should be computed not down to 
Ri, (which would introduce much faster oscillations) but down to ~ r/j -|- Ar^. 

The H-wiggles have not been observed in the probabilities for other objects. In the 
case of SMA solution of the solar neutrino problem, due to smallness of the vacuum mixing 
angle the wiggles are so small that they are simply unobservable. For the LOW solution 
the mixing at the production point is very strongly suppresed. For supernova neutrinos, 
again the amplitude of the wiggles is very strongly suppresses and, moreover, the period is 
so small that the wiggles are averaged out. 

3.8 L-wiggles 

The L-wiggles are the consequence of mixing of vim and V2m in a given fiavor state Va and 
the interference of the transition amplitudes 

Va Vim Vi and Va V2m ^i, (3.35) 



where i = 1,2. For energies above the H-resonance the states z^im and 1^2™ are mixed mainly 
in and I'r, and their admixture in is neghgible. Therefore the L-wiggles are realized in 
channels with initial states of and u^-. (see Fig. |3|, upper panel). The transitions ( 3.35| ) 
appear due to adiabaticity violation in the L-resonance. 

For definiteness we will describe the wiggles in the channels in the case of normal 
mass hierarchy. According to ( p. 17 ) in the production point has the following matter 
eigenstate content 

(3.36) 







Note that below the H-resonance energy the admixture of j/2m is suppressed by small- 
113. In contrast, for E > Ej^ we have 8^3 ~ 1 and according to ( [3.36 ) 



ness of S12, 



-C23l^lr 



5236 l^2r 



Consider evolution of this state. Its I'lm and z^2m components undergo various transforma- 
tions with the following amplitudes in factorization approximation: 



\A{iyir 

\A{iyir 



V2) 

^1) 



Pl , 



\A{u2r, 
\A{v2r. 



T^2m 
V2m 



^2 



)| = V(i-^'h)(i-Pl) 



vx)\ = V(l - Ph)Pl . 



Furthermore, at the resonance crossing A{h'im i^2m) = —A{i'2m i^im)- (For com- 
pleteness: |A(f2m — ^ 1^3)1 = VPh, which is irrelevant for this interference.) Note that vim 
crosses only L-resonance, whereas i'2m crosses both resonances which is indicated by two 
arrows. Summing up the amplitudes from different channels of transitions we obtain 



P*{i'/i 

P*{i^^i 
P*ii^ti 



^1) 

2^2) 
J^s) 



C23Vr^ - e-''t'^S23e'Wi'^-PH)PL 

C23 + e-'^^S23e'Wi^-PH){l-PL) 
sir? 923PH ■ 



(3.37) 



Here (j)L is the phase difference of the amplitudes collected from the production point to 
the end of L-resonance region where the adiabaticity is restored. To a good approximation 
one can use for the adiabatic phase difference: 



dx {H2m - Hi 
r 



(3.38) 



Ro 



The probabilities ( 3.37| ) can be rewritten as 



P*{i^^i 
P*{v^i 
P*{1^^l 



yi) = 4(1 - Pl) + sUl - Ph)Pl - I, 
^2) = cI,Pl + sl,{l - Ph){1 - Pl) + I, 
^3) = 43PH, 



(3.39) 



where the interference term equals 



= sin 2^23 cosi^L + <5)7(1 - Ph)Pl{1 - Pl) 



(3.40) 



Note that the CP-violation phase enters together with the oscillation ( "strong" phase) and 
to disentangle the former one needs to know (p^. Let us first summarize properties of the 
interference term. 

1. The interference term is not suppressed by small mixing in contrast to the H- wiggles 
considered in the previous section. Here the admixture is determined by the vacuum mixing 
agle 023 which is close to maximal. In fact, in the energy range where Pl ~ 1/2, this term 
can have the amplitude of the order 1. 

2. The interference term appears in the transitions to vi and 1^2 and not to u^. It has an 
opposite sign in — vi and — )■ 1^2 probabilities (as it should be according to unitarity) . 

3. The interference term vanishes if Ph = 1 which, in turn, is realized for ^13 — t- 0. Non- 
zero 1-3 mixing is a necessary condition for its appearance. The amplitude is maximal 
if the transition in H-resonance is adiabatic. With increase of energy the adiabaticity in 
H-resonance is broken, Ph — )• 1, and this suppresses the H-wiggles. 

4. A necessary condition for appearance of wiggles are P^ ^ and Pl ^ 1, that is, the 
adiabaticity should be broken in the L-resonance but it should not be broken very strongly. 

5. The energy region where the L-wiggles are realized (more precisely, the lower border of 
this region) is determined by the following two conditions: 

• the mixing of the eigenstates uim and i'2m should be large enough which happens in 
the region of the H-resonance and above it; 

• the adiabaticity in the L-resonance should be broken. For the density profiles we are 
discussing the adiabaticity starts to be broken in the energy range of H-resonance 
and above it (which is to some extent accidental). As a result, the L-wiggles appear 
in the same energy region as the H-wiggles. 



The flavor probabilities can be obtained plugging expressions ( 3.39 ) in (|3.6|) . In par- 
ticular, 



P{U^ ^ l^e) = COS2 023 [\Uel\^ " PUPell^ " |f/e2|')] + 



Sin2 023(1 - Ph) [\Ue2\^ + PLi\Uel\^ " |C/e2p)] + (|f/e2p " \Uel\^)I, 



Let us consider separately the interference terms in the flavor probabilities. According 

to i ^M ) 



P{V^ Ve)int = {\Ue2? - Wel?)!^,^ 
Pij^f, I^m)™* = (|t^M2p - \Uf,l\'^)If„ 
Pil^f. ^ l^rUt = i\Ur2\^ - \Url\^)I^. 

For the selected values of the mixing ang les = 0.5 and sin^ 26*13 = 0.08 we have 
P{i^^ ^ lyeUt = -0.372/^ , P(l/^ ^ UrUt = 0.318/^ , PiUf, ^ U^),nt = 0.054 



in agreement with the results of numerical computations (see Fig. |3|, lower panel). The 
wiggles in the — >• z^e and — >• Vr channels are large and in opposite phase and their 



amplidude decreases with increase of energy because P/, — )■ 1 and Ph — >• 1- The L- wiggles 
are suppressed in the — >■ channel. 

Let us now consider the phase of oscillatory behavior ( 3.38| ). The phase is collected 



(i) from the production point rj to the H-resonance, r|j . Here the difference of eigenvalues 
is essentially given by the vacuum frequency: Am1-^/2E; (ii) in the H-resonance region 
and below: here the difference of eigenvalues start to decrease down to ~ Am^2/2-£'j (iii) 
in the L-resonance, where the frequency is ~ Amf2/(2-E sin 2^12). The phases collected in 
regions (ii) and (iii) are small. In the L-resonance the oscillation length becomes larger 
than whole baseline. For instance for E = 200 GeV, we obtain 1^ ^ Q ■ 10^^ cm ^ rj^. 
Although now the phase is collected from much larger region than in the case of H-wiggles: 
from Tj to + Arj^ (e.g. 2- 10^^ cm, for E ~ 200 GeV), the phase (j)L is comparable to (j)H 
(collected by 0.7 • 10^^ cm). The reason is that the phase between + Arjj and + Ar^ 
is relatively small. Consequently, the period of L-wiggles is comparable to the period of 
H-wiggles but certain phase shift is present. 

3.9 Dependence on the density profile 

The dependence of the probabilities on characteristics of density profile of the star (the 
initial density no (see Fig. ^) and gradient of density change in the envelope) is shown in 
Fig. ^ for profiles A and B. 

Let us first consider dependence of probabilities on k. With decrease of k the adia- 
baticity determined by A„ = r^/k becomes stronger. Indeed, for fixed energy the resonance 
layer has larger radius r/j as k decreases. Consequently Ph becomes smaller and 

• the 1-3 dip (in v^, — )• fg channel) becomes deeper and reaches the adiabatic minimum, 
Pmin = sfs) even for small values of 1-3 mixing (see Fig. |5|, upper panel); 

• adiabaticity is broken at higher energies {E > 170 GeV for sin^ 26I13 = 0.08), and 
correspondingly, the wiggles (which are manifestations of the adiabaticity breaking) 
shift to higher energies; 

• the amplitude of L-wiggles does not change significantly; 



• period of wiggles decreases. This follows immediately from Eq. ( p. 34 ): is larger 
and therefore the phase for a given E increases. 

Similar features are present for — )• u/^ modes (see Fig. ^, lower panel). With decrease 
of k the adiabaticity both in H- and L- resonances improves. This means that in the range 
E > E^ the probabilities reach their adiabatic values. The amplitude of wiggles becomes 
smaller and period is substantially smaller. The latter is due to increase of rj^ as in the 
case of H-wiggles. 

For antineutrinos, decreases of gradient leads to stronger matter effect in the range 
(lO'^ — 10^) GeV, essentially extending region of the adiabatic conversion in the L-resonance 
to higher energies. 

Effect of the inner density increase is illustrated in Fig. |3|. There are two consequences 
of this increase: (i) shift of the resonance energies and therefore whole energy region of 



matter effects to smaller energies; (ii) improvement of adiabaticity in both resonances due 
to the fact that for a given energy the radii of resonance layers increase. This, in turn, 
leads to the fohowing observational consequences: 

• the dip due to 1-3 resonance shifts to smaller energies and starts at = 60 — 80 GeV 
for the factor of 2 density increase; 

• the dip (peak) reaches the adiabatic minimum (maximum) for the survival (transi- 
tion) channels; e.g. for P(fe — )■ i^e) ~ sin^ and P{Ufj, — )• I'e) ~ 0.5(1 — sin^ Oi^); 

• the H-wiggles have smaller period and so there are more wiggles on the nonadiabatic 
edge; 

• the L-wiggles have smaller amplitude at low energy which then increases with energy. 
This is related to good adiabaticity in the L-channel and small Pl which increase 
with energy. Period of wiggles becomes smaller. 

• Due to better adiabaticity the asymptotics are achieved at higher energies. 

The antineutrino probabilities are affected very weakly and mainly in the range (10^ — 
10^) GeV. 

On the contrary, with decrease of the initial density in the envelope and increase of 
the gradient, the matter- affected region shrinks: the region shifts to higher energies but its 
non-adiabatic edge - to lower energies. 



4. Properties of conversion probabilities 

Let us consider dependence of conversion probabilities in specific channels on neutrino 
parameters ^13, 623 as well as on the type of mass hierarchy. 

4.1 Probabilities of u^, — up transitions 



In Fig. |4.l| we show the Ue ^ u/s, (3 = e, ij,,t probabilities as functions of neutrino energy 
for different values of the 1-3 mixing angle. Properties of these probabilities can be well 
understood using the considerations in the previous sections. 

1. In asymptotics according to ( 3.13 ) we have the averaged Su— oscillation probability in 
vacuum: 

P(l/e ^Ue) = J2Pei\'' = COS^ ^13 (l " ^ 2^12^ + sin^ ^13. (4.1) 

For z/e — )• Ufj, channel the probability equals P{i'e — ^ i^n) = |f^ej|^|f^/xjp) etc. 

2. In the range Ej^ <^ E <^ E^, inserting the matrix of probabilities ( p.l9[ ) into Eq. ( 3.1C| ) 
we obtain 

P(l/e ^ Ue) = \Ue2\\l " |t/e3|') + Peil^ = |C/e2|' " | f/e3 1 ' ( I f/e2 P " |C/e3p), (4.2) 

or in terms of the angles: i-*(fe — >■ i'e) = cos^ ^13 sin^ ^12 -|-sin^ ^13. For the — ?• channel 
we have 

m ^ I^m) = \U^.2\Hl - \Ue3f) + |^7^3p|t/e3p = \U^2\^ + | ^7e3 P (| C/^3 1 ' " \U^2f), (4.3) 
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Figure 5: The — > vfi (top panel) and — > vfi (bottom panel) transition probabilities as 
functions of the neutrino energy for two different stellar density profiles, p cx r^^ (red curves) and 
p oc r~^^/^ (black curves), but for the same initial density at Rq. We used sin'^ 26*13 = 0.08, Sep — 
and normal mass hierarchy. 



and similar expression for P{i^e — ^ i^t) with substitition — s- iC/rip. 

3. For E > and large enough 1-3 mixing the adiabatic evolution gives according to 

(|3!20|) and (|3lc| ) 

P{iye ^ l^e) = \Ue3\^ P{l^e ^ '^^^) = \U^3\^ P(l/e ^ Z^r ) = | f/r3 1 ' • (4.4) 

These results reproduce the probabilities at ~ 200 GeV for the 1-3 mixing sin^ 20i3 > 
0.1. For smaller 1-3 mixings the adiabaticity is broken already at -E ~ E^ and therefore 

If neutrinos are produced not too far (in energy scale) from the H-resonance, the 
fg— state contains non-negligible admixture of the z^2m eigenstate: 

fe = cos 6'°3 Z^2m + sin 6*13 I'Srn ■ (4.5) 

Adiabatic evolution of this combination will give then 



P^Ue 1^1) = 0, P.iUe V2) = COS^ e\^, P,{ue ^ 1^3) = sin^ 



13) 



and instead of ( |4.4D for i/^ — )• i/^ probability we obtain 

P{i^e i^e) = cos^ cos^ 6*13 sin^ 6*12 + sin^ 6'?3 sin^ 6*13 . 
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Figure 6: Probabilities for Ve — > (top panel) and — > (bottom panel) transitions for 
different values of 1-3 mixing: sin^ 20i3 = 0.04 (black lines), 0.08 (red) and 0.15 (blue lines). We 
take profile A, Sep — 0, sin^ 023 = 0.5 and normal mass hierarchy. 

For E > with increase of energy the adiabaticity violation becomes important. 

The analytic results presented here allow one to understand the dependence of the 
probabilities on ^13 and 623. In asymptotics and in the intermediate region (plateau) the 



probabilities only weakly depend on ^13 (see Fig. iA_): the corresponding corrections are 



proportional to \Ue3\ = sf^ (see (p~l|), ([4.2[), ( [4. 3D ). The strongest effect is in the range of 
H-resonance and above it, where the 1-3 mixing is enhanced. For the dip in the adiabatic 
case we have P 

^ (|4.4D . With decrease of s^g the adiabaticity becomes broken and the 



survival probability increases. In the limit of very small 5^3 the second dip disappears and 
we will have only one dip due to 1-2 mixing. The probabilities in other channels have an 
"inverted" dependence. 

The dependence of probabilities on the 2-3 mixing is shown in Fig. ^. Note that 
the Ue — >• fe probability does not depend on 623- The dependences of other probabilities 
have rather interesting feature: there are two energies Ei ~ 100 GeV and E2 ~ 800 GeV 
at which the probabilities do not depend on ^23 7 and moreover, at E2 ~ 800 GeV the 
probabilities in all the channels are equal. With increase of S23 the Ue i^fj. probability 
decreases at E < Ei and E > E2 and it increases in the interval Ei — £'2- Indeed, 



according to 
(0): Pii^e - 



in the intermediate range P(fe — f/^) = |C^/i2p ~ C23, and in the 1-3 dip 
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sL- In asymptotics, a weaker change can be immediately 
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Figure 7: The Ve — > (top panel) and v^t_ — > vp (bottom panel) probabilities as functions of 
energy for different values of 2-3 mixing: sin^ 023 = 0.4 (blak curves), 0.5 (red curves) and 0.6 (blue 
curves). We used profile A, sin^ 26*13 — 0.08, 6cp — and normal mass hierarchy. 



understood from expression P(z/e — )• ^ |C^eiP|f^^ti |^ + 1 142 P I t^/i2 P • The strongest 
dependence is in the plateau. The dependence for — )• Vt channel is just inverted. 

4.2 Probabilities of i/^ up and Ur — 5- J^/? transitions 

The probabilities — )• up) as functions of neutrino energy for different channels and 
two different values of 1-3 mixing are shown in Fig. 4.1. The main difference from the 
Ug— channels is that is not the eigenstate of the propagation basis (in contrast to 
Ue) but combination of these states. This leads to more complicated expressions for the 
probabilities and additional interference effects. The probabilities equal 

i 

and explicit expressions for probabilities in different energy regions are given below. 
The asymptotic values of probabilities for E <^ El and E » El equal 

P{u^^up) = Y,\U^.i\^\Upif■ (4.6) 

i 

In the plateau we obtain from (3.10) and ( |3.19| ) 



or 

\2\TT |2 \TT |2|rr |2|rr |2 



PC. . ., ^ _ l^rsl \Uf)l\ , I^msI I^e3| 1^/32^ , ut \2\tt |2 ^4 
^(^M ^ - i_|f;^3|2 + T^^^^^p + 1^/^31 |f/.3| , (4.7) 



where we used that 



sm 623 = TTTlo^ ^23 = r. rrrri- (4. 



Above the H-resonance for large 1-3 mixing, which satisfies the adiabaticity condition, we 
find from (|3lo| ) and (|3^) 

(|f//32|VM3|' + |f^/3l|Vr3|'). (4.9) 



1 - |f/e3P 



It reproduces in the first approximation correct values of the probabilities in the dip at 
E ~ (150 — 170) GeV for large 1-3 mixing and for small gradients. Certain deviation of 
numerical results from analytic ones is a manifestation of the adiabaticity violation. For 
maximal 2-3 mixing we obtain from ( [4.91) : P(f^ — s- vp) = 0.5(1 — |C/^3p). 

As in the case of Ve — ^ the change of 1-3 mixing mainly affects the probabilities 
in the energy range above the H-resonance. With decrease of ^13 the jump probability 
increases and according to ( 3.40| ) the amplitude of wiggles, oc ^Jl — Ph, decreases. 

The v-T-— amplitudes and probabilities (they appear at the border of envelope due to 
oscillations inside jets) can be found from the z/^— amplitudes obtained in the previous 
subsection by substitutions: — )• C/t-3i sin 6*23 — ?• cos ^23, cos023 — ^ — sin023. 

4.3 Probabilities in antineutrino channels 

There is no resonances in antineutrino channel for normal mass hierarchy and the depen- 
dence of the probabilities on energy is simpler (see Fig. ^) . The 1-3 mixing is not enhanced 
and the conversion effects are mainly due to large 1-2 mixing at low energies (where this 
mixing is not suppressed). 

The asymptotic values of the probabilities are the same as in the neutrino channels 
(|4.6| ). There is no CP-violation asymmetries. In the intermediate region {E^ < E < E^) 



we obtain from (|3.10D and (|3.22| ) 

m ^ = |C^/3ll'(l - |C^e3|') + \UfS3\^\Ue3\^ = \U(,i\^ - | C/e3 1 ' (I 1 ' " \Up3\^), 

P{yf. ^ i^/3) = Y3^7^ {\Ur3\^Wp2\^ + \U^,3?\Ue3?\Upi\^) + | ^7^3 1 V/33 1 ' , (4.10) 

(/3 = e, /i,r). They reproduce well the result of numerical computations in the range 
(10 - 50) GeV (see Fig. |). 

In the range E > E^, (150-180 ) GeV, the 1-3 mixing is suppressed by matter, so 
that ~ 1, s^^ ~ in ( [3.22[ ), and the flavor transition probabilities in the adiabatic 
approximation are 



/31l 



|2 



P{y, ^ = (If^r3p|?7^2p + I^^3P|?7^3|') • (4.11) 



Comparing the first probability with the one in ( [4.10| ) we conclude that P{i^e ~^ ^e) only 
slightly increases in comparison with the intermediate region. With further increase of 
energy the probability decreases due to the adiabaticity violation, approaching the vacuum 
oscillation result. Correspondingly, P(Pe — ^ ^t) is "inverted". 

The probabilities in antineutrino channels (in contrast to the neutrino resonance case) 
rather weakly depend on energy. Their dependence on the 1-3 mixing is weak (Fig. |8|): 
according to ([4.10| ), the corrections to P{i'e — ^ ^13) are of the order jC/eSp, and corrections 
to P{v^ — )• Up) are even smaller. The probabilities differ from their asymptotic values 
mainly in the intermediate range. Interestingly, the probabilities with in final state (at 
least for zero 5 and maximal 2-3 mixing) are not affected by matter. 
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Figure 8: The same as in Fig. 4.1 for the antineutrino channels. 



According to Fig. ^ the — )• fe probability does not depend on 6*23. Dependences 



of probabilities on 2-3 mixing are not very strong and well described by Eqs. ( 4.10 ) and 



4.4 Probablities for the inverted mass hierarchy 

Now the L-resonance is in the neutrino channels, whereas the H-resonance is in the an- 
tineutrino channels. Therefore the probabilities are given by some combinations of the 
results obtained for v and v for the normal mass hierarchy (see Fig. |l^). The differences 
appear when interplay of the two resonances (in the case of NH) becomes important. 
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Figure 9: The same as in Fig. ^for the antineutrino channels. 

Let us consider first the neutrino probabilities, (i) The asymptotic values are the same 
as in the NH case, (ii) In the intermediate region we can use the same result as for NH 
( ]4.3| ) and ( [4.2[ ). (iii) Above the H-resonance according to and ( p^ we have 

|2 



P{iy, ^ up) = \U, 



/32| 



(/3 = e, ^l, t), 



(4.12) 



which should be compared with the results in Eq. (4^) for the normal mass hierarchy. 
In particular, P(fe — ^ z^e) = |f^e2p- Comparing this result with (4.4) we see that the 
probability slightly increases in the H-resonance in contrast to the decrease in the NH case. 
According to Fig. 12 for sin^ 2^13 = 0.08 the adiabaticity is broken already at E 
therefore values of probabilities deviate from those in Eq. ( [4.12 ). 



and 



In the intermediate energy region for the (f^ — )• u^) channels we have the same results 
as in the case of normal mass hierarchy ( 14.7] ). Above the H-resonance, again the difference 
from the NH case appears since there is no level crossing. Using (3.23) we obtain flavor 
probabilities 



'23 



l-|t/e3|' 



) 



which should be compared with Eq. ([4.S|). 

Let us consider the antineutrino channels. In the intermediate energy range for the 
initial Ug we have the same results as in the NH case ( 4.10| ). Above the H-resonance (which 
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is now in the antineutrino channel) using Eq. ( 3.24| ) we obtain 



P{y^^yp) = \Up:i\\ (4.13) 
1 

1-We3| 



They reproduce correct values of the probabilities in the interval ii^ ~ (150 — 170) GeV. 
4.5 Effects of oscillations inside jets 

In Fig. |l^ we compare the probabilities with and without oscillations in jet. As follows from 
these plots, the jet effect is very small and appears in the region above the H-resonance. 
This can be explained in the following way. Oscillations between neutrino production 
point inside jet and inner part of the envelope are described by Aj^tiya t'^) introduced 
in Eq. (^^). We assume that in the neutrino production region a jet has an average density 



1.5 • 10 cm . The corresponding resonance energies equal: E'^- = W GeV and 
E^- = 6-10^ GeV for the 1-2 and 1-3 mass splits correspondingly. Total matter width 
(column density) equals njTy ~ 9 • lO'^^ cm~'^ which is much smaller than the refraction 
length and therefore matter effect can be neglected. Indeed, for E < Eji the potential 
is smaller than the kinetic (vacuum) term, whereas for E > Eji the oscillation length is 



much smaller than the baseline, so that vacuum mimicking situation is realized [31|. In 
general, amplitude of vacuum oscillations in the flavor basis can be written as 



AjetK ^ = + UUU'{s [e'^'^'' - l) + U'^2U'{2 [e'''^'' - l) , (4.14) 

where </>3i = nx/l^i and 021 = '7rx//2i are half oscillation phases inside jet. Vacuum 
oscillation lengths for the two modes equal /13 = lO^'^ cm (£'/100 GeV), and /12 = 3 • 
10^^ cm (E'/lOO GeV). The typical length of the neutrino trajectory inside jet r^, ~ 6 • 
109 cm, and therefore for £" > 60 GeV the phase 021 ~ 0. So, we can neglect oscillations due 
to 1-2 mass split and consider non-averaged vacuum oscillations due to 1-3 mass splitting 



only. Then the amplitudes ( 4.14 ) become 



where 



ri{x)= e^^-^^i^^^-l . (4.16) 
According to ( |3.4D the mixing matrix elements in the propagation basis are given by 

Thus, the S-matrix of transitions inside jet ( [4.15 ) can be written as 



1 + r/sfg r?si3Ci3 
Sjet = I + vUsUl = I 1 



77S13C13 1 + r/cfg 
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Figure 10: The — > I'p (top panel) and — > (bottom panel) transition probabilities with 
oscillations in the jet (thick red curves) and without oscillations in jet (thin black curves). We used 
sin^ 2013 — 0-08, sin^ 6*23 — 0.5 and Scp = in normal mass hierarchy. 



Inserting this matrix into ( |3.11| ) we obtain the total S-matrix inside the star: 



( ' 



+ VS23 

13Cl3S23e' 



??Sl3Cl3C23e 



r/si3Ci3 

C23 S23e'\l + V43) 
-S23 C23e 



X 5, 



(4.17) 



''(1+^4) 



where iSenv)ai = A^nvil^'a ^ '^i)- 

In the adiabatic approximation the S-matrix in the envelope equals the mixing matrix 
at the inner part of the envelope: Senv = ^'^d the latter is given in Eq. (3.16). Inserting 



this matrix into ( 4.17 ) we obtain explicitly 



C?3(l + r]sl^) - ??S?3Si3Ci3 S?3(l + r/S^g) + r?C?3Si3Ci3 

-C23 f?C?3S 13 0135236*"^ " S?3S23e*'^(l + ^/C^g) ?7Sl3'Sl3Cl3'S23e*^ + C?3S23e*'^(l + rjcj^) 

-S23 VCi3Si3Ci3C23e'^ - s?3C23e*'^(l + rjc^^) r/s?3Si3Ci3C23e*'^ + c?3C23e"^(l + 770^3) , 

(4.18) 

Consider this matrix in specific energy ranges. In the intermediate energy range (C13 = C13, 



'13 



S13) it is 




Cl3 
-Sl3'S23e 
-Sl3C23e 



•513(1 + ??) 
Cl3S23e'\l + v) 
Cl3C23e'\l + v) 



Corrections due to oscillation inside jet are given by the terms oc rj. Since (1 + 77) = e'^"^^^ 
and 1(1 + 7y)p = 1, no corrections to the probabilities appear due to oscillations in jet. An 
absence of the corrections is related to the particular initial state (mixing in matter) in the 
envelope, the adiabatic evolution in the envelope and loss of coherence on the way from 
the star to the Earth. 

Above the H-resonance {ci^ ~ 0, s^s ~ 1) have in the adiabatic approximation 

/ -??si3Ci3 1 + r/sfg \ 

S,= i-C23 -5236^^(1 + r?cf3) r/Si3Ci3S23e'M • (4.19) 
\-S23 -C23e'Hl + Ticjr^) ??Sl3Cl3C23e'V 

The probabilities are given by moduli squared of the 5* elements: P^.{va ^i) = \iS*)ai\'^- 
Then, according to ( ^.191 ) corrections appear. Note that ll+rycfgp = \ l+r]*sl^\^. Explicitly 
we obtain from ( |4.19D the probabilities for fg— channels: 

P*il^e l^l) = 0, P*il^e 1^2) = S^cLl^P' P*i''e ^ i^s) = |1 + S?3^P- (4-20) 

Now corrections are non-zero but they are suppressed by small factor sf^. 
For the f^— channels we obtain from ( 4.19 ) 

P*{^t^ J^i) = C23, P*{i^f, U2) = S23II + c?3^P> P*{^^l ^3) = S23■Sl3C?3kl^ (4-21) 

and again the corrections are small, being suppressed by 8^3. 

From expression for S^: it follows that in the adiabatic case the probabilities do not 
depend on S: the S^:— matrix elements either do not depend on 6 or are proportional to 
the overal phase factor e^^ . 

Similarly, for antineutrinos the probabilities of transitions are not affected by oscilla- 
tions inside jet in the intermediate range. Above the H-resonance, according to ( 4.19D , 

P^{i?e ^ 1^1) = \l + sj^-ql"^ , P*(z^e j^2) = 0, P*(z^e 1^3) = SigCiglr/l^. 
P*(j^At Z^l) = S?3Cl3S23l^l^ ^'*(l^M ^2) = C23, (P^ -)> 1^3) = ^^3 1 1 + cf gT/^ , 

(compare with the corresponding results for neutrinos). As for the neutrino case, here we 
have non-zero corrections which are suppressed by small factor sf^. 

In the case of inverted mass hierarchy the amplitudes of oscillation inside jet are the 
same as for normal mass hierarchy, at least in our approximation (2z^— vacuum oscillations). 
So, the amplitudes Aj^t equal the amplitudes in Eq. ( [4.15 ). Consequently, general formulas 
for the total amplitudes inside a star are given in Eq. ( [4.17 ). 

For intermediate energy range we have exactly the same results as in the case of NH: 
there is no corrections due to oscillations inside jet. The difference appears at energies 
above H-resonance since now the level crossing and mixing in initial state are changed (see 
Fig. We obtain the following expressions for probabilities: 

P*(z^e l^l) = 0, P*(z^e ^ 1^2) = \1 + sfgr?!^, P^{Ue U3) = Sl3Cl3kl^' 
P*(j^/. ^ Z^l) = C23, P*(^'m 1^2) = S23S?3Ci3|t/|^ P* (l/^ ^ 1^3) = S23 | 1 + CigT/^ . 



Note that probabilities of transitions to i^i are the same as in NH case, whereas the prob- 
abihties for transitions to 1^2 and ^3 have been interchanged. The conclusion is the same 
as for NH: corrections are non-zero but they are suppressed by small factor sfg. In the 
antineutrino channels in the intermediate energy range results are the same as in the case 
of NH. Above the H-resonance all the probabilities coincide with those we had obtained in 
Eqs. ( [1.20 ) and (4.21) for neutrinos. 



An overall conclusion is that there is no corrections due to oscillations inside jets in 
the energy range between resonances and small corrections appear above the H-resonance. 
In general, the jet effect appears if the adiabaticity is broken. 

4.6 CP-violation effects 

Let us consider dependence of the oscillation probabilities on the CP-violation phase (see 
Figs. |ll] and |l^). Outside the matter affected energy interval, E = {1 — 10^) GeV, the 
probabilities equal the averaged oscillation probabilities in vacuum. Important feature 
is that the coherence of the mass states is lost after evolution inside star. Effects of the 
{va — ^ J^i) transitions sum up incoherently. So, projection back to final flavor state depends 
on moduli and no interference effects appear. These probabilities depend on moduli 

squared of the mixing matrix elements and the interference terms are absent. In this case 
no CP-odd effects appear, the CP- as well as T-asymmetries vanish and P(fQ, — t- i^^) = 
P{ua ^p)- The dependence on CP-violation phase originates from \Ucii\'^ {a = ^jl,t, and 
i = 1, 2) and is rather weak since the terms with 5 are proportional to sin0i3. 

The interference and CP-violation can be related to projection of the initial state onto 
the propagation basis states. To describe this let us first neglect oscillations inside jet and 
consider specific conversion channels. 

1. Vf> —7- Ufi channels: since v^, is the component of the propagation basis, no CP- 
violation appears in the projection and further convertions. P(z^e — ^ i^e) does not 
depend on 5 (see Figs. |ll| and |l^ top panels). P^Ue — >• Vji) (/3 = ^, ''") do depend on 6 
via |C//3ip (i = 1,2) (see Eq. ( p.27D ) via projection of the mass states onto the flavor 
state in the detector. Furthermore, for maximal 2-3 mixing with change of 5 the 
(i^e i^fi) probability is transformed into {ve — ^ ^t) and vice versa, when 5 increases 
from to vr, and aX 6 = 7r/2 

2. — )• channels: now according to ( |3.3D 

= C2^v^ + 5236*^1^;, (4.22) 
and the transition z/^ — )• i/j has two channels via v'^ and z/^ which can interfere. So, 



as follows from ( [1.22| ) the probability equals 



2 



PA^i, Vi) = C2zA^{y'^ Ui) + S23e'M*(i/; ^ Ui) . (4.23) 

The interference term depends on 6. As we saw before in the adiabatic case A^,{u'^ — )• 
Vi) = 5ii and A^{v'^ Ui) = (0, -s?3, c^g) (see Eq. (|3T^)). 
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Figure 11: The probabilities as functions of energies for different values of the CP-phase Sep — 
(black curves), 7r/4 (red curves), and tt/2 (blue curves). We take normal mass hierarchy and 
sin^ 2^13 = 0.08. 

Let us consider the interference in various energy ranges. In the intermediate energy 
range, according to ( p.l6| ) u'^ ~ —uim and v'^ ^ —Si3i'2m + cisi^sm- Since in this range 
z^im evolves adiabatically to i^i and the latter is orthogonal to the rest of the state, no 
interference appear between v'^— and z^^— channels in ( 4.23 ), and consequently, 



I 2 I 



This is in accordance with results ( ^.181) . So, the interference and dependence on S require 
adiabaticity violation, which is realized above the H-resonance. For E > we have 



'13 



1, f' PS —I'lm, and u'^ ^ — f2m- Therefore 



and due to adiabaticity violation, both amplitudes in this equation are non zero simulta- 
neously for i = 1 and 2. The phase 6 appears in the interference term in combination with 
the oscillation phase and therefore with change of 5 the wiggles above the H-resonance 



shift (see Fig. 11). It might be practically impossible to disentangle the effect of 5 from 
(f)L- Similar consideration holds for Vr- 

Thus, in adiabatic energy range the dependence of the flavor probabilities on 5 follows 
from the projection of mass states onto the flavor states in the detector only. The projection 
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Figure 12: The same as in Fig. ^ for antineutrinos and the case of inverted mass hieararchy. 



of the original flavor state onto the eigenstate of propagation gives an additional dependence 
on 6 only if the adiabaticity is broken inside the star. In this connection let us take into 



account oscillations inside jets. From Eq. (4.17) we have 



A^{l^e l^i) = Aenvil^e t'i) + f?Sl3 [sisAenvit^e t^i) + CrsAenviK ^i)] , (4-24) 

where the second term is the jet effect. No dependence on the CP-phase appears here. For 
the — )• Ui transitions, according to ( 4.17] ), 



+ S23Cl3e'^r] [si3Aenv{j^e ^^i) + Ci3Aenv{K ^i)] , (4-25) 



where the same combination of the amplitudes as in (|4.24D enters the correction term. 
For rj = (no oscillations inside jet) this expression is reduced to the one in Eq. (|4.23| ). 
The dependence on the CP-phase here is explicit. All other factors and amplitudes do not 
depend on 5. 

Let us consider the antineutrino channels. General form of the amplitudes is the same 
as for neutrinos with the following changes: (5 — >• —6, and 0™ — )■ 6^ which corresponds 
to change V — >• —V. The amplitudes of oscillations inside jet are unchanged; they are 
essentially the 2^— oscillations in vacuum in our approach. 



5. Fluxes at the Earth 



5.1 Flavor ratios 

We compute the fluxes at the Earth separately for neutrinos and antineutrinos keeping 
in mind that in future experiments the signs of charged leptons may, in principle, be 
determined. The fluxes at the Earth equal: 

where q = e, r, and e is given by Eq. (^). Similar expressions will be used for antineu- 
trinos with substitution e — >■ e. Let us introduce the flavor ratios as 

^ _ P{v^l yg) + ePjVe 

"■"/^ ~ ~ P{y^ ^ u^) + eP{ve ^ ' 

and for antineutrinos 

normalizing all the fluxes at the Earth to the f^— flux at the Earth. 

According to Fig. |l] the ratio of original fluxes decrease with energy as 

^1^, r 1 ; < 10 GeV 



$0 1 0.5 ; ^> 300 GeV. ^^'^^ 

The ratios contain complete information relevant for observations. Properties of the flavor 
ratios can be easily undertstood from properties of probabilities. If the charge of lepton 
produced by a neutrino is not defined, and thus, the experiment sums up signals of neutrinos 
and antineutrinos, observables are determined by the ratio: 

where ^q, describe ratios of antineutrino and neutrino cross-sections and corresponding 
efficiencies of detection. Note that contribution of antineutrinos is suppressed by about 0.2 
- 0.3 due to smaller cross-section (^q, ~ 1/2) and smaller original flux (another factor 1/2). 

At energies below 10^ GeV, and Vr can not be distinguished: both produce the 
showering events. The double-bang events which are signatures of the i^,-— interaction can 



be identified at much higher energies |34]. Thus, with existing detectors one can study the 
ratio of the showering and tracking events (with muon track in the final state) r^/j/^^. In 
the first approximation it is determined by the ratio of the sum of fg— and Vr— fluxes to 
the f^— fluxes: 

This expression should be corrected (even neglecting misidentification of events): (i) certain 
part of z/'t-— flux contributes via transitions: z/,- — > t — )• to the traking events (see e.g. [^]); 



(ii) neutral current interactions (not affected by oscillations) contribute to the showering 
events. 

To a good approximation, for E > 100 GeV, the charged current cross-sections are 
equal for all three flavors and therefore we can sum up the and Uj- fluxes at the detector. 
(Note that still one should subtract the channel Vr ^ t ^ which leads to the tracking 
event.) Therefore if the initial flux is composed of i/^, the total flux which produces 
showering events is determined by the probability 



Psh = ^ ^^){\Ue^\^ + (1 " h^,)\U, 



Tl I ) 



where ~ 0.2 is the fraction of r — t- decays. Using unitarity we find 

Psh = ^P*(l^^ ^ Vi){l - |C/^^|2) = 1 - ^n(z^^ ^ ^^^)|f/M^|' ^ 



2 

Tl I ) 



and the latter equals Psh = ^ — Ptoti^iJ. i^fj) — b^ P*{vpL — ^ ^^jJl^riP- If dominates in 
the initial state, Ptotii^fj. t-'^i) determines the probabilities of all observable events when 
r — )• /X decays are neglected. In the presence of the original z^e— flux we have 

Psh^l- Pii^f. ^ + e[l - P{ve ^ y^,)]- 
5.2 Properties of the flavor ratios 

In Figs. 13 - |l^ we show the flavor ratios r^,/^ and r^/^ (a = e,/i, r) as functions of 
neutrino energies, values of neutrino parameters and density profiles. For illustration we 
use two extreme original flavor contents: (i) e : 1 : 0, where e = e{E) and according to 
Fig. II it decreases from 0.5 at £^ < 10 GeV down to 10~^ at high energies. This content 
would corespond to stong suppression of muon decays in medium of jet without secondary 
acceleration, (ii) 1:2:0 which corresponds to e = 0.5 and does not depend on energy. 
This could be realized when additional acceleration of secondary muons occurs in jets. 
Consider first e : 1 : 0. For very small e we have 

Furthermore, since Pii'fj, — >• i^f^) has weak dependence on energy, the ratio essentially 
repeats behavior of -P(f^ — )• i^a)- According to Fig. |l3| (upper panel), and in agreement 
with our analytic considerations, the ratio r^/^ has an asymptotic value ~ 0.7 and the 
deviation of this value from the flavor equilibration is due to the original flavor content. In 
the intermediate region (plateau), due to equality Pi^i'^ — )• fe) ~ Pi^fi — ^ i^fj,) for maximal 
2-3 mixing we have r^/^ w 1. The 1-3 peak (at E = 200 GeV) is slightly enhanced in 
comparison to the peak in probability. 

In the plateau for maximal 2-3 mixing we have P^i'^ — ?• i't) < Pi'^ii ~^ T^fi)^ and conse- 
quently, r^i^ ~ 0.8. The wiggles of P{i'^ — )• v^-) and r^j^ have an opposit phase as compared 
with wiggles of P{i'^ — )• fe) and ^g/^, and the amplitude of wiggles of showers/tracks ratio 
is suppressed. For maximal 2-3 mixing the asymptotic value equals r^/^ ~ 1. Due to 
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Figure 13: Neutrino and antineutrino flavor ratios (fluxes normalized to the z-'^ flux) as functions 
of the neutrino energy for two different values of 1-3 mixing, sin^ 26i3 = 0.04 (blue lines) and 0.08 
(red lines) and two different original flavor contents: e : 1 : (upper panel) and 1:2:0 (lower 
panel). The thick (thin) lines correspond to the neutrino (antineutrino) fluxes. We used profile A, 
sin^ 023 = 0.5, Sep — and normal mass hierarchy. 



opposit phase the r^/^ and r^/^ wiggles compensate each other in the ratio r^j^/i^., leading 
to substantially smaller amplitude of Vgh/tr wiggles (see Fig. p^). 
For anineutrinos (see Fig. 13) the ratio 



follows to a large extent the ratio of original fluxes in ( ^.1[ ) in the case of NH. An additional 
distortion due to oscillations is rather small. In particular, difference of high and low energy 
asymptotics in flavor ratios is related to difference of ratio $p^/$|^^ at high and low energies. 

For the original flavor content (1:2:0) (bottom panels) the fg— flux and its trans- 
formations substantially change the flavor ratios at the Earth. In the intermediate range, 
still Tg/^ w 1 as a consequence of the equality P(fe — )• Ve) ~ Piye ^ii) for maximal 2-3 
mixing. At and above the energy of H-resonance the probability P{ye — ^ ^e) has a dip, 
whereas P{ye f/x) has a peak and above the peak/dip they have wiggles of the opposit 
sign. Therefore the wiggles of ratio rg/^ are strongly attenuated. This ratio deviates from 
1 (equilibration) which is the signature of the matter effect only in the region (10^ — 10^) 
GeV and maximal effect is about 30%. Similarly, for the ratio r^/^ substantial attenuation 
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Figure 14: The same as in Fig. |l^ for inverted mass hierarchy and only for sin^ 26is — 0.08. 



occurs. Moreover, the wiggles of r^/^ and r^/^ have opposite phase, and therefore the 
wiggles in the Tg/j/^^ are further suppressed (see Fig. p^). 

As follows from Fig. 13, the highest sensitivity to ^13 is in the range of H-resonance 
peak and wiggles: E = (10^ — 10^) GeV. (Recall that in the f^— channels the peak is also 
due to resonance enhancement of the 1-3 mixing.) With decrease of ^13 adiabaticity is 
broken stronger and the amplitude of wiggles becomes smaller. 

Let us consider the case of inverted mass hierarchy (see Fig. 14). For the (e : 1 : 0) 
original content, as in the case of NH the asymptotics is r^/^ ~ 0.7 and in the plateau 
^e/^j ~ 1- In the H-resonance region there is no peak, the ratio decreases with increase of 
energy from 1 to the asymptotic value. This dependence is modulated by small wiggles. 
The peak appears in the H-resonance region in the antineutrino ratio r^/^. 

For the original flavor ratio (1:2:0), due to strong compensations of contributions 
from different channels the matter effect on the flavor ratios is small. 

Let us consider dependence of the flavor ratios on other neutrino parameters. 

1. The phase of wiggles changes with 6: wiggles shift (see Fig. 15). Changes in plateau as 
well as in the asymptotics are rather weak. 

2. Dependence of flavor ratios, as functions of energy, on the 2-3 mixing is strong (see 
Fig. |l^. The strongest dependence is for r^/^ (NH) in the plateau region: as we already 
mentioned, for maximal 2-3 mixing the probabilities of (i'^ — t- fe) and (f^ — t- f^) transitions 
coincide. However with change of 623 they change in the opposite way: one increases and 
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Figure 15: Neutrino and antincutrino flavor ratios as functions of the neutrino energy for two 
different values of the CP phase S: (red lines) and 7r/4 (blue lines) and two different original 
flavor contents: e : 1 : (upper panel) and 1:2:0 (lower panel). The thick (thin) lines correspond 
to the neutrino (antineutrino) fluxes. We used profile A, sin^ 6*23 — 0.5, sin^ 813 — 0.08 and normal 
mass hierarchy. 

another decreases, as can be seen in Fig. |l^. The — t- u-j- probability changes much 
weaker. Note that asymptotic values also vary rather substantially, whereas the change in 
the 1-3 peak is small. Thus, determination of r^/^ in the plateau region would be the most 
sensitive to search for deviation of the 2-3 mixing from maximal. For the original flavor 
content (1:2:0) the dependence on 6*23 is weaker. 

Consider dependence of the flavor ratios on the density profile. With decrease of 
density gradient, k, (see Fig. |l^) the evolution becomes more adiabatic, the 1-3 peak shifts 
to smaller energies and becomes wider; the amplitude of wiggles decreases, the region of 
wiggles extends to higher energies. 

For the original flavor content (1:2:0), there is substantial cancellation of contri- 



butions from the z^„ 



Uf. and I'p 



transitions, since the transition and the survival 



probabilities have opposite dependence on energy. 

Dependence of the ratios on the initial density no is strong (see Fig. 18): with increase 
of riQ the resonance peak shifts to lower energies, the period of wiggles decreases and 
number of wiggles increases; the amplitudes of wiggles become slightly smaller. In the case 
of (1 : 2 : 0) original content, there is strong cancellation of contributions of (t'^ — t- v^) 
and {ue — )• i^e) especially in the region of 1-3 peak. As a consequence, the region of wiggles 
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Figure 16: Neutrino and antineutrino flavor ratios as functions of the neutrino energy for different 
values of sin^ 6*23: 0.4 (black lines), 0.5 (red lines), 0.6 (blue lines), and two different original flavor 
contents: e : 1 : (upper panel) and 1:2:0 (lower panel). The thick (thin) lines correspond to 
the neutrino (antineutrino) fluxes. We used profile A, sin^ 2^13 = 0.08, Sep = and normal mass 
hierarchy. 



starts from about 200 GeV. 

In general, in the case of original (e : 1 : 0) content, dependences of the flavor ratios 
on the neutrino parameter as well as on parameters of the star are much stronger, and this 
makes hidden jets to be more prospective than the observable jets as far as investigation 
of neutrino parameters is concerned. 

In Fig. |l^ we show the ratio of shower-to-track events on the neutrino energy defined 
in Eq. ( |5.2| ). Quahtatively, this dependence repeats the corresponding dependence of the 
flavor ratio r^/^ with L- wiggles and asymmetric asymptotics at high and low energies. 
However relative size of the wiggles, and in general, matter effects (deviations from the 
averaged oscillation result) is further suppressed. As we mentioned before, this suppression 
is due to contribution of the Vr transition, r^/^, as well as the contributions of 

antineutrinos which have small matter effect for the normal mass hierarchy. According 
to Fig. 19, the structures in r^/i/j^ due to conversion in the matter of the star are up to 
40 — 50%. This should be compared with factor of (2 - 2.5) effect in the ratio r^j^. The 
effect is smaller for the inverted mass hierarchy. Substantial deviations from the averaged 
VO result are in the energy range (102 -2-103) GeV which can be extended to (30 - 5 • 10^) 
GeV for smaller density gradient or larger density in the inner part of the envelope (bottom 
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Figure 17: Neutrino and antineutrino flavor ratios as functions of tire neutrino energy for different 
density profiles: A (red lines) B (blue lines) and two different original flavor contents: e : 1 : (upper 
panel) and 1:2:0 (lower panel). The thick (thin) lines correspond to the neutrino (antineutrino) 
fluxes. We used sin^ 29i3 = 0.08, Sep = and normal mass hierarchy. 

panel). The size of the matter effect is much smaller: about ±5% in the case of initial ratio 
flavor ratio (1:2:0). 

6. Discussion and conclusion 

Hidden jets realize unique set up (configuration and conditions) in which at one hand 
protons can be accelerated up to high energies and then produce high energy neutrinos, 
and at the other hand these neutrinos cross large enough column density of matter on the 
way from production region so that flavor conversion in the matter of the star is important. 
We present detailed and comprehensive study of flavor conversion of neutrinos from hidden 
sources (jets) - astrophysical sources invisible in non-thermal 7 rays or X rays. Our results 
differ from the results of paper [|l6|. 

1. The conversion is affected by matter of a star and the Earth in the energy interval 
~ (1 — 10^) GeV. In this range the probabilities and flavor ratios substantially deviate 
from probabilities and ratios given by the averaged vacuum oscillations. The borders of 
this interval are determined by the largest density in the envelope. El = £'i(no), and 
gradient of density in the envelope ~ i?^Am|^/47r rather than on R-^, the radius of the 
stellar envelope. The interval expands for larger initial density and smaller gradient. For 
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Figure 18: Neutrino and antincutrino flavor ratios as functions of the neutrino energy for the profile 
A with different inner number densities of electrons no = 10^"^ cm~^ (red lines) and uq — 2 ■ 10^^ 
cm~^ (blue lines) and two different original flavor contents: e : 1 : (upper panel) and 1:2:0 
(lower panel). The thick (thin) lines correspond to the neutrino (antineutrino) fluxes. We used 
sin'^ 2013 — 0.08, Sep = and normal mass hierarchy. 



energies below this interval matter effects can be neglected and with decrease of energy the 
probability converges to the averaged VO probability. Above the interval the adiabaticity 
is strongly broken and flavor evolution inside a star can be neglected. With increase of 
energy the probability converges to the averaged VO probability again. 
2. We discused in details physics involved in the production of neutrinos. Neutrinos are 
produced in vr— , /i— and iiT— decays in strong magnetic fields which leads to very short wave 
packets of neutrinos. This, in turn, results in quick separation of the wave packets and loss 
of coherence, in spite of high neutrino energy. The coherence length is comparable with 
radius of star. The coherence is not restored in a detector due to finite energy resolution 
and large baselines. Therefore the overal evolution consists of the flavor-to-mass states 
transitions inside the star, loss of coherence, propagation of the nmmass staes in vacuum; 
so that incoherent fluxes of mass states arrive at the surface of the Earth. Detector projects 
these mass states back onto flavor states incoherently. 

Inelastic scattering becomes important for high energies ~ 10^ GeV. To a good approx- 
imation this scattering is flavor invariant, and therefore factors out and does not influence 
the flavor conversion picture. It can be included in determination of the original neutrino 
spectra. We have estimated that some effect of — scattering may show up at high 
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Figure 19: Ratio of the shower-to-track events as the functions of the neutrino energy for two 
different initial flavor content: e : 1 : (blue lines) and 1:2:0 (red lines) upper panel: for normal 
(sohd) and inverted (dashed); we used profile A; lower panel: for tiq = 2 x 10^-^ cm^^ (sohd lines) 
and for smaller gradient k = 17/7 (dashed lines). We used sin^ 26'i3 — 0.08, Sep — 0. 

energies where usual flavor conversion effect becomes small. 

3. Mainly the neutrino flavor change is due to adiabatic and partially adiabatic conversion 
in the envelope of a star. Oscillations inside the Earth lead to an additional distortion 
of probabilities below 10 GeV. We show that probabilities as functions of neutrino energy 
have several generic features which are well controlled by the density profile crossed by 
neutrinos (initial density, density gradient). This includes 

• Plateau in the intermediate energy range which is due to adiabatic conversion in the 
1-2 resonance region. The Earth matter effect produces an oscillatory dip in the 
plateau at < 10 GeV. 

• The dips or peaks (depending on channel) related to the 1-3 resonance at -E > Eh- 

• Non-adiabatic edge of the energy profile is modulated by wiggles. The wiggles are new 
dynamical feature which is not realized in other objects. The wiggles are manifesta- 
tions of interference induced by adiabaticity violation in the H-resonance (H-wiggles) 
and in L-resonance (L-wiggles). 

Different features of the energy profile of the matter effect such as low energy bor- 
der, plateau, the 1-3 peak/dip (its position and size), wiggles and asymptotics depend on 



parameters of neutrino and star differently, so that one will be able to disentangle these 
dependences in principle. Position of these features depends on density profile of the star: 
in particular, on the density at the border between jet and envelope (no) and gradient of 
density in an envelope (k). The size of the matter effects depends on values of neutrino 
parameters as well as on original flavor content. In the adiabatic range there is no depen- 
dence on the profile (apart from lower border of the regions). Therefore plateau as well as 
the 1-3 peak (for large enough 6*13 or small gradient) do not depend on a profile. Above 
the H-resonance in the adiabaticity violation range (the wiggles region) characteristics of 
probabilities depend both on profile and on neutrino parameters. In particular, number of 
wiggles is a measure of the density gradient in the envelope. 

4. At the Earth one can, in principle, determine the flavor ratios ^e/^) '^t/^i '^e/fj.j ^■f/A" 
^p.//! as functions of neutrino energy. To a large extent these ratios reproduce the energy 
dependence of probabilities. Independent measurements of these ratios would give very rich 
information on neutrino parameters, density profile of star and original fiavor ratios which 
encode information about conditions inside jets. Effect of neutrino flavor conversion on 
these ratios is strong, modifying the ratios in certain energy ranges by factor 2-3. With 
present experimental techniques and detectors, however, one can determine the ratio of 
shower-to-track numbers of events only, without separation of neutrinos and antineutrinos. 
Sum of contributions from neutrinos and antineutrios as well as from i^g and Ur (which 
have an opposite phase) lead to significant damping of the observable conversion effect in 
fsh/tr- The effect depends strongly on the original flavor content. It can reach 30 - 40 % 

for the content e : 1 : (e ^ 1) and it is about 5 10% for 1:2:0. In reality one can 

expect some intermediate situation. 

5. Observability of the neutrino signal from hidden jets is beyond the scope of this paper, 
and here we just add some comments. There are two issues related to observability: (i) 
number of event in a given detector, and (ii) a possibility to identify signal and extract 
it from a background produced by the atmospheric neutrinos. These issue have been dis- 
cussed partly in |l^ . Detection of signal from individual source would require large detector 
and/or a source in a nearby galaxy such as M82 and NGC253. Signatures include high en- 
ergy neutrino events from certain directions during ~ 10 sec. which can be repeated. There 
should be a correlation between the jetted neutrino signal {E >1 GeV) with burst of low 
energy thermal neutrinos (ii^ ~ 5 - 50 MeV) generated during the core collapse. Thermal 
10-50 MeV neutrinos are also produced in a jet since the shell density and temperature 
are initially very high at the base of the jet. No significant delay is expected between the 
thermal and high-energy neutrinos that are produced in the jet. There can be significant 
delay between the emission of core-collapse thermal neutrinos and jet thermal/high-energy 
neutrinos depending on whether the core collapses directly to a black hole or via an inter- 
mediate neutron star stage. The thermal neutrino emission from core-collapse preceding 
the high-energy (and more thermal) neutrino emission and the time-delay between the two 
can be used as a probe to learn about the core-collapse and jet formation processes. 

In case of a successful supernova, with or without jet break out, optical lightcurve can 
be extrapolated back to the time of explosion and to check temporal correlation with any 
neutrino signal detected from the same direction in the sky. This can largely reduce the 



atmospheric neutrino and muon background. A neutrino detector setup outside the Earth's 
atmosphere, e.g., in space or on the Moon, as the technology develops, can of course greatly 
improve the detection prospect of astrophysical neutrinos. 

It is expected that parameters of jets change during jet duration which leads to time 
dependence of neutrino signals we have discussed. In principal, it opens unique possibility 
to monitor evolution of jets with neutrinos. 

The neutrino oscillation signatures from hidden jets are modified in the case of isotropic 
diffuse flux. To identify this flux one can use distortion of the energy spectrum in the energy 
range where atmospheric neutrinos are not affected by oscillations. In general if the diffuse 
flux is dominated by hidden sources as discussed here, the measured flavor ratio will differ 
significantly from the usual (1:1:1) ratio expected from the optically thin sources in the 
GeV - TeV range. 
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A. Meson and lepton energy losses 

High energy vr and K are subject to energy losses due to hadronic irp and Kp interactions 
with cross-sections (T-^p/Kp ^ (2/3)o"pp, and due to synchrotron radiation in the strong 
magnetic field before they decay. The cooling time scales due to these processes are 

for hadronic interactions, where AE'^^j^^ ~ 0.5i?^^^^ is the average inelasticity and we 
ignored the logarithmic factor in cross-section, and 

/ ^/ N -1 / N -2 ( 0.34 s ; j = K 

^ ^ ^ ^ t 7.0 X 10--^ s; j = fi . 

for the synchrotron cooling. As can be seen, from these equations, the meson energy losses 
are initially dominated by hadronic and later by electromagnetic process |, |. With 
a decay time scale t'j = Tj{Ej/nij), where tj is the mean lifetime, n and K decay 
without significant energy losses below E'^^ ~ 400 GeV and E'^ ~ 50 GeV (which 
corresponds to t'^^^ ^ ^'had)- "^^^ electromagnetic losses become important at energies 
above £^^^2 ~ "^^^ ^'n b2 ~ ^^0 GeV (which corresponds to t'^^ t'^^^). Thus 
we can define a suppression factor to be multiplied with the production fluxes of vr and K 
before they decay as 

V , , , 1 ^-^^i 

Cj{Ej) ^ <j t'j,had/^'j,dec ^ (-^j/^j-.w)-^ ; ^j,bl ~ ^ - ^j,b2 



Note that muons from tt- or X-decays are also subject to synchrotron energy losses which 
dominate above 5 ~ 6 GeV and the //—flux before decay will be suppressed by a factor 
~ {E'^/E'^^)^'^. Thus, the z^'s from /i— decay channels contribute little to the total flux 
at high energy. Also, the contributions to the fluxes from np— and i^Tp— interactions is 
small and we ignore those. 



B. Comparison of the conversion formulas in Ref. ||16|] and in this paper 



General formula ( |3.6D for the conversion probabilities differ from that used in the paper 
fl^ . Consequently, results of this paper differ from the results in [16|. According to [16| 
the flavor probability equals 



(B.l) 



where v-y is a flavor neutrino state, and Pv{T^'y — ^ ^p) is the averaged oscillation probability 
in vacuum. Taking explicit expression for the latter we can rewrite the probability ( |B.1| ) 
as 



or 



I "-^ 7« I ) 



(B.3) 



where A* 



is the amplitude of probability of the corresponding flavor transition. 



On the other hand, expression for the probability 
as 



we use in this paper can be rewritten 



^7 ) ^"yi 



(B.4) 



which clearly differs from (B^). The difference is that instead of the flavor states at the 
exit from the star, we take the mass states which are the eigenstates of propagation in 
vacuum and these states lose coherence. In ( |B.3| ) the averaging was taken for oscillations 
outside the star, and oscillations of mass states are not averaged inside the star. Therefore 
Eq. ( |B.3D leads to fast oscillatory picture which corresponds to oscillations inside the enve- 



lope in the interval from rjet — R*- In our consideration the wiggles also appear. However 
they appear in the adiabatisity violation range only, i.e. above the H-resonance. There is 
no wiggles in the adiabatic part. Furthermore, the period of wiggles in the energy scale is 
much larger: the phase of wiggles is determined by the oscillation phase on the way from 
the production point to the H- resonance region (for H-wiggles) and to the L-resonance 
(for the L-wiggles) which is much smaller than total size of the envelope rjet — -R*- 



C. The Earth matter effect 

The mass-to-flavor transition probabilities inside the Earth can be written as 

Pe{i^^ ^ 1^0) = m23rsSEU')pi\^ = \ {U23rsSEUi3Ui2)pi\\ (C.l) 



where according to (|3.4| ) U' gives projection of the mass states onto the states of the 
propagation basis, Se is the S-matrix for transitions inside the Earth in the propagation 
basis, and U23rs projects the propagation basis onto the flavor basis. (If Se = I, the 
probabiUty (CJ) is reduced to \{UpMNs)i3i\'^ ■) 

Let us find exphcit expression for Se for E > 1 GeV where substantial neutrino flux 
from jets exists. For the Earth densities this energy range is far above the L-resonance and 
the 1-2 mixing is strongly suppressed. So, the problem is reduced to 2z/— oscillations due 
to the 1-3 mixing. For definiteness we will take normal mass hierarchy and consider for 
simplicity the evolution in the approximation of constant density. In this case 



(C.2) 



where 

D = diag{e-'^'',l,e-''^'^). (C.3) 

Here (I)i2 = 4>\ — (t>2-, 4'32 = 4'3 ~ 4'2, and = Hit are the oscillation phases of the 
eigenstates in matter i^im (Hi are the eigenvalues in matter). The mixing matrix in matter 
in the propagation basis, U^, for E ^> El equals, according to (3.5), 
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Using this expression for and 
Earth: 
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the evolution matrix in the 



(C.4) 



A, 



m2 I m2 -4032 
^13 ~r '^13 ^ 



(C.5) 



Inserting ( |C.4| ) into ( C.l ) we obtain the matrix of flavor-to-mass transition probabilities. 
In particular, we have 



Then for the flavor probability f, 
and (|1|): 



Si3Ci2^eTp, 



lee — •Sl3Si2^ er I ) 



J^e) = \ci3C12Ae 
I'e) = \ci3Si2Ae 
Ve) = \s13Aee + Ci2,Aer\'^ . (C.6) 

; — )• i>e with all transitions included we obtain using 



Cl3Sl2|ci3^ee " Sl3^eTp + s{^\si3Aee + Ci^Ae 



(C.7) 



This is reduced to the probability without oscillations inside the Earth ( |4.2| ) for A^ 
and Aer 



0. In the first approximation in S13 we have from (C.7) 



''12l^ee| 



■512 



1 - sin^ 2df^ sin^ 



923 



(C.8) 



According to this result the plateau is modulated by resonance oscillations: the oscil- 
latory dip appears in the range E ~ E^^ ~ 6 GeV. Typical width of the dip is given by 
AE ~ 2£'(3 tan 2^13. So, for sin 2^13 = 0.08 the modulations are in the range (4 — 8) GeV. 
Maximal depth depends on the zenith angle of the neutrino trajectory. The result ( |C.8 ) 
differs from the Earth matter effect on the solar neutrinos where both plateau and the 
Earth effect are due to 1-2 mixing and the Earth matter effect enhances the survival prob- 
ability. Here the effect is opposite: the probability becomes smaller. The difference stems 
from the differentce of density profile of the Sun and an envelope of a star. In particular, 
density in the central part of the Sun, where neutrinos are produced are much higher than 
maximal density in an envelope. As a result, in the case of jetted neutrinos, the Earth 
matter effect due to 1-2 mixing is at very low energies (< 0.2 GeV) - essentially in the low 
asymptotic region - below E^. The Earth matter effect in the 1-2 mixing plateau is due to 
1-3 mixing. This (to some extend) is similar to the 1-3 dip at E > 10^ GeV. 

Detection of the Earth matter effect could give an information about neutrino prop- 
erties and also about direction to the star. It is not clear though, if this effect can be ever 
extracted from the atmospheric neutrino background. 
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